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An exact expression is derived for the electron Green function in two-channel Kondo models with 
one and two impurities, describing the crossover from non-Fermi liquid (NFL) behavior at interme- 
diate temperatures to standard Fermi liquid (FL) physics at low temperatures. Symmetry-breaking 
perturbations generically present in experiment ensure the standard low-energy FL description, but 
the full crossover is wholly characteristic of the unstable NFL state. Distinctive conductance line- 
shapes in quantum dot devices should result. We exploit a connection between this crossover and 
one occurring in a classical boundary Ising model to calculate real-space electron densities at finite 
temperature. The single universal finite-temperature Green function is then extracted by inverting 
the integral transformation relating these Friedel oscillations to the t matrix. Excellent agreement 
is demonstrated between exact results and full numerical renormalization group calculations. 



PACS numbers: 75.20.Hr, 71.10.Hf, 73.21. La, 73.63.Kv 



I. INTRODUCTION AND PHYSICAL PICTURE 

The full power of the renormalization group (RG) con- 
cept is perhaps most clearly seen in its application to 
quantum impurity systems.^ The classic paradigm is the 
Kondo model, ^ being the simplest to capture the fun- 
damental physics associated with all quantum impurity 
models: universal RG flow from an unstable fixed point 
(FP) to a stable one on successive reduction of the tem- 
perature or energy scale. The Kondo model describes a 
single local spin-^ 'impurity', coupled by antiferromag- 
netic exchange to a single channel of noninteracting con- 
duction electrons. Here, perturbative scaling arguments"^ 
indicate an RG flow from a high-energy unstable 'free 
fermion' FP (describing a free impurity decoupled from 
a free conduction band), to a low-energy stable 'strong 
coupling' FP (where the impurity is screened by conduc- 
tion electrons via formation of a 'Kondo singlet'). This 
RG flow is characterized by a scaling invariant — the 
Kondo temperature Tk — which sets the crossover en- 
ergy scale. But analysis of the crossover itself goes be- 
yond simple scaling ideas and the conventional RG pic- 
ture. Wilson's numerical renormalization group'' (NRG) 
allows an exact nonperturbative calculation of certain 
thermodynamic and dynamical quantities which show the 
crossover (for a review, see Ref. 5). Universal scaling of 
all physical quantities in terms of the crossover scale Tk 
confirms the basic RG structure of the problem. 

However, a different RG flow occurs when the im- 
purity is coupled to two or more independent conduc- 
tion channels.^ In this multichannel Kondo model, the 
frustration inherent when several channels compete to 
screen the impurity spin renders the strong coupling FP 
unstable. A third FP at intermediate coupling^ then 
dictates the low-energy physics. This FP exhibits non- 
Fermi liquid (NFL) behavior, including notably a residual 
entropy^ of ifcs ln(2) in the two-channel Kondo (2CK) 
model. The crossover from the free fermion FP to the 
2CK FP has been the focus of much theoretical atten- 
tion. In particular, solution of the model using the Bethe 



ansatz yields the exact crossover behavior of thermody- 
namic quantities;^ while NRG has been used to calculate 
thermodynamics® and dynamics^'^" numerically. It was 
also shown recently that this 2CK physics can arise in 
odd-membered quantum dot rings^"' and chains, ^^ and in 
quantum box systems. ^"^"^^ 

Indeed, the same type of NFL behavior^^"^^ arises in 
the two-impurity Kondo (2IK) model. ^° The tendency 
to form a trivial local singlet state is favored by an ex- 
change coupling acting directly between the impurities; 
while the coupling of each impurity to its own conduction 
channel favors separate single-channel Kondo screening. 
The resulting competition gives rise to a critical point^*^ 
that is identical to that of the 2CK model with additional 
potential scattering. ^^ 2IK physics is also expected to ap- 
pear in certain double quantum dot systems, ■^^ and other 
even-impurity chains.^* 

A description of the NFL FPs of such two-channel 
models in terms of an effective boundary conformal field 
theory (CFT) shows that the operator controlling the 
FP has an anomalous scaling dimension. ^^'^'^ This im- 
plies unconventional energy/temperature dependences of 
physical quantities such as conductance Gc{V,T), mea- 
sured as a function of bias voltage V and temperature T. 
In the 2CK device of Ref. 24, ^/V/T^ and ^/tJT^ cor- 
rections to the NFL FP conductance predicted by CFT 
were directly observed in experiment. Similar signatures 
are expected^®'^^ in the channel-asymmetric 2IK model; 
although leading linear behavior emerges in the symmet- 
ric 2IK.^^ This behavior is of course in marked contrast 
to {VlTKf and {T/TkY Fermi liquid (FL) behavior ob- 
tained ubiquitously in the single-channel case.^^ 

The NFL FP itself (and the crossover to it) has now 
been rather well studied. However, NFL physics is ex- 
tremely delicate: various symmetry-breaking perturba- 
tions destabilize the NFL FP^^'^'^ and generate a new 
crossover scale T* . At T = 0, the impurities are thus 
completely screened and all residual entropy is quenched. 
Indeed, regular FL behavior,^ including the standard 
{V/T*Y and (T/T*)"^ corrections to conductance, must 
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Figure 1: Schematic phase diagram for tlie 2CK and 2IK 
models, as a function of temperature T, external energy scale 
w, and symmetry-breaking perturbation strength S. The three 
FPs of each model give rise to three distinct regimes: free 
fermion, NFL and FL. We considered the NFL to FL crossover 
at T = in Ref. 26, indicated by arrow (a). Here we generalize 
the results to finite temperature, arrow (b). 



appear at low temperatures T <C T* and energies V <Si 
T*. Therefore, no evidence of nascent NFL physics can 
manifest in the immediate vicinity of the FL FP. Only on 
fine-tuning the perturbation strength (5 — > to the crit- 
ical point so that T* -^ does one obtain NFL physics 
on the lowest energy scales. 

But RG analysis in the vicinity of the free fermion, 
NFL and FL FPs implies two successive crossovers, with 
Tk setting the energy scale for flow to the NFL FP, and 
T* characterizing flow away from it. Even in the FL 
phase away from the critical point (which is the generic 
case relevant to experiment), NFL behavior can be ob- 
served at higher temperatures and energies, provided 
there is good scale separation T* <^ Tk (see Fig. 1 for 
a schematic phase diagram). Li this case, conductance 
Gc{V, T) through the 2CK quantum dot device of Ref. 24, 
or proposed 2IK devices, ^^'^^ should exhibit a clean NFL 
to FL crossover. 

In Ref. 26, we considered this conductance crossover 
at T = as a function of bias V , corresponding to the 
crossover labelled by arrow (a) in Fig. 1 (and V here 
playing the role of the external energy scale ui). It was 
shown^^ that the full crossover is wholly characteristic 
of the high-energy NFL state. The T = crossover is 
expected to describe the behavior at very low temper- 
atures. From a scaling perspective, RG flow is cut off 
on the energy scale F - T, so Gc{V, T) ~ Gc{V, 0) for 
T <^T* since there are no further crossovers below T* . 

By contrast, at higher temperatures T :^ T* (<C Tk), 
no evidence of the NFL to FL crossover will be observed 
(see Fig. 1), and only the NFL FP is probed. Indeed, this 
is the likely scenario in the experiment of Ref. 24: rather 
than tuning to the critical point 6 — 0, signatures of the 
true FL ground state are simply washed out by tempera- 
ture. But the behavior as a full function of temperature 
T and energy scale w for finite perturbation strength S is 
much more subtle, and naturally strengthens connection 
to experiment. Exploring the third temperature axis in 



Fig. 1 , and considering the resultant NFL to FL crossover 
[eg. arrow (b)], is thus the focus of the present work. 

In this paper we combine Abelian bosonization 
methods^'^'^'^'^^ with the powerful machinery of CFT^^'^^ 
to obtain an exact description of the NFL to FL crossover 
in two-channel Kondo models. In particular, we calculate 
the full electron Green function at finite temperature, 
from which conductance follows. ^^''^° The field theoretic 
description links the 2IK model with a classical Ising 
model on a semi-infinite plane. ^'^ Application of a bound- 
ary magnetic field h in this boundary Ising model (BIM) 
results in RG flow from an unstable FP with free bound- 
ary condition /i = to a stable FP with fixed boundary 
condition h — > ±00.^^''^^ This RG flow due to h is iden- 
tical to that occurring between NFL and FL FPs in the 
2IK model due to a small perturbation 6.^^ Indeed, an 
emergent symmetry of the NFL FP in the 2IK model, ^^ 
together with the common CFT description of 2IK and 
2CK models, ^^'^^'^"^ implies the existence of a single uni- 
versal NFL to FL crossover function for both models, re- 
sulting from any combination of relevant perturbations.^^ 

Exact results'^^'^"' for the BIM are the source of our 
solution, which becomes exact when there is good scale 
separation T* <C Tk, as sought experimentally. 

The exact crossover Green function at T = was calcu- 
lated in Ref. 26 by exploiting the above connection. How- 
ever, ambiguities appear at finite temperature which pre- 
vent straightforward generalization of those results. Thus 
we take a different route here: the BIM solution is used to 
calculate real-space Friedel oscillations around the impu- 
rities at finite temperature, which are themselves related 
by integral transformation^^ to the Green function. The 
problematic analytic continuation is avoided in this way. 

The paper is organized as follows. In Sec. II we intro- 
duce the 2CK and 2IK models, together with representa- 
tive symmetry-breaking perturbations that generate the 
NFL to FL crossover. We then present and discuss our 
main results for the exact finite-temperature Green func- 
tion along the crossover. The corresponding conductance 
crossover for quantum dot systems that might realize 
2CK or 2IK physics is then calculated. In Sees. Ill- VI 
we derive the analytic results. First we consider the 2IK 
model at T = with a single detuning perturbation. In 
Sec. Ill we calculate the resulting crossover Green func- 
tion, exploiting the analogy to the BIM. In Sec. IV we 
extend the calculation to finite temperatures, extracting 
the desired t matrix from Friedel oscillations. The results 
are generalized to the 2CK model in Sec. V and to an 
arbitrary combination of perturbations in Sec. VI. Ex- 
act results are compared with finite-temperature NRG 
calculations in Sec. VII. Other quantities showing the 
crossover such as entropy and nonequilibrium transport 
are then briefly considered in Sec. VIII. The paper con- 
cludes with a general discussion in Sec. IX, and details 
of certain calculations can be found in the appendices. 



II. MODELS AND RESULTS 

We consider the standard 2CK and 2IK models, 

H2CK = Hq+JS ■ (soL + Sob) + SH2CK, (1) 

H21K = Ho+J{Sl ■ SQL + Sr ■ Son) + KSl ■ Sji + 6H21K, 

(2) 

where Ho = J2a,k''k^kcra'^kaa describes two free con- 
duction electron channels a = L/R, with spin density 

SOa = Eaa'^'Ljl^^<y')%a'a (^^nd ^La == EfcV'La) 

coupled to one spin- 2 impurity S (2CK) or two impurity 
spins Sl,r (2IK). For SH2CK = 0, the NFL ground state 
of H2CK is described by the 2CK FP. Likewise, a critical 
inter-impurity coupling Kc can be found such that the 
ground state of H21K is again a NFL,^" and is similarly 
described by the 2CK FP for 5H2ik = O.^^-i^ 

Relevant perturbations to the above models (embodied 
by 6H2CK and 5H2ik) are those that destabilize the NFL 
FP, and result in a FL ground state. A new scale T* 
is thus generated, characterizing RG flow from NFL to 
FL FPs. The relevance of such perturbations can be 
traced to the breaking of certain symmetries, ^^'^"^ such 
as parity or particle-hole symmetries. In fact, there are 
many possible perturbations to the 2CK and 2IK models; 
but two perturbations may be considered 'equivalent' if 
they break the same underlying symmetry — and their 
effect on the low-energy physics will be identical. ^^'^^ 

For concreteness, we consider now the simplest pertur- 
bations which exemplify such symmetry-breaking, and 
which in combination generate all possible NFL to FL 
crossovers at low energies/temperatures. Specifically, 
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(3) 

describes L/R channel asymmetry in the 2CK model for 
A^ ^ 0, while charge transfer between the leads is em- 
bodied in the A^; and Ay components of the first term 
[here r (a) are the Pauli matrices in the channel (spin) 
sector]. The second term describes a magnetic field act- 
ing on the impurity. For the 2IK model, the critical point 
is destabilized by finite (Kc — K), and also through 
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Figure 2: Schematic illustration of possible 2CK (a) and 2IK 
(b) setups to measure conductance. The left lead is 'split' into 
source and drain, allowing the resulting conductance through 
the attached impurity to be measured. 



themselves be different'^^). Thus we do not consider such 
perturbations explicitly here. 

A. Quantities of interest 

Signatures of the NFL to FL crossover on the scale 
of T* should appear in all physical quantities. In 2CK 
or 2IK quantum dot devices which could access this 
physics, ^^'-"^^'^^ the quantity of interest is the dl/dV 
conductance G';^{V,T) = {2e^h-^G^)G'^iV,T) through 
channel a = L ov R. Here, C^ = 4r^rjJ/(r^ + V^)^ 
describes the relative strength of coupling to source and 
drain leads (see Fig. 2 for an illustration of the setup). At 
zero-bias V = Q, the conductance is given exactly by,^^ 

dw t„a[^,T), (5) 



G?(o,r) = 



^E 



duj 



1] ^ is the Fermi function, and 



where f{uj,T) = \e}^l^ 

taa (w, T) is the energy- resolved local density of states (or 

'spectrum'). 



t„c,{y>,T) = -T\v Im7^Q_^Q(a;,r), 



(6) 



itself related to the t matrix"'^ 7^0^0-0,(0;, T), describing 
scattering of a cr =t or J, electron within channel a = L 
or _R, with bare lead density of states per spin v. Note 
that for equal hybridization to source and drain leads, 
F" = FJJ, Gq = 1 is maximal; while in the asymmetric 
limit F^ < rj^, G;J < f . In the latter case, the weakly- 
coupled source lead probes the system perturbatively, so 
the system remains near equilibrium, even at finite bias 
1^ > 0. The resulting conductance is then simply,^*^ 



iil2IK = ^(FLfl,VLL^O<xi? + H.C.) + B.-(S'L-S'fl,), (4) G^^{V^T) ^hY, f 



duj- 



-df{Lu-V,T) 
duj 



t„o.{uj,T), (7) 



where Vlr describes electron tunneling between the leads 
and Bs the application of a staggered magnetic field. 
Channel anisotropy could also be included in the 2IK 
model, but the critical point can always be recovered-*^^ 
on retuning K. Similarly, spin-assisted tunneling be- 
tween channels Yla '^Oah'^OaR^L ■ Sr + H.c. (as arises in 
a two-impurity Anderson modeP^) is expected to have 
the same destabilizing effect as the Vlr term in Eq. 4, 
since they both have the same symmetry at the NFL 
FP'^^ (although the resulting crossover energy scales may 



where t^ra is the equilibrium (zero-bias) spectrum. 

The t matrix itself must show signatures of the NFL 
to FL crossover since scattering is purely inelastic at 
the NFL FP,^^'^'^ but inelastic scattering must cease at 
energies <^ T* , where the impurity degrees of freedom 
are fully quenched.^* Thus the crossover also shows up 
in conductance. Our goal here is to calculate the full 
crossover t matrix, and hence conductance, at finite tem- 
perature for the 2CK and 2IK models in the presence of 
symmetry-breaking perturbations described by Eqs. 1-4. 



B. Survey of results and discussion 



with Xf = {X2, A3, Ai}. For the 2IK model it is. 



In the next sections we derive an exact expression for 
the desired t matrix, describing the universal crossover 
from NFL to FL behavior in the 2CK and the 2IK models 
at finite temperature — and which as such generalize the 
results of our previous work in Ref. 26. Here we pre-empt 
the full derivation, and present our key results. 

The NFL to FL crossover is characterized by a low- 
energy scale T* arising due to the presence of symmetry- 
breaking perturbations to the 2CK and 2IK models. It 
is given generically by^^ 



'<2IK 



T* = A^ 



(8) 



where A^ = X],=i -^?- The eight contributions correspond 
to relevant perturbations which have distinct symmetry 
at the NFL FP. Two perturbations which have the same 
symmetry correspond to the same Xj . The perturbations 
given in Eqs. 3 and 4 are classified viz, 



Table I: Classification of perturbations 



A. 


2CK model 


2IK model 


Ai 


ciuA.VTk 


c^{K,~K)/VT^ 


A2 


ciuA^Vr^ 


cvRe i^VlrVTr 


A3 


ciuAyy/Tii 


cylm uVlrVTk 


As 


cbB/VTk 


cbBs/VTk 



where Xb = {A5,A^,A^} = {A4,A5,A6}. The pertur- 
bations associated with coupling constants A7 and Ag do 
not conserve total charge, -^^'^^ and so are ignored here 
(although we note that such operators can be of impor- 
tance, for example, in the context of strongly correlated 
superconductors'"' ) . 

ci,cv,cb — C(l) are fitting parameters'^^ which de- 
pend on the model and on J. and Tk cx e~'^ is the 
Kondo temperature, characterizing RG flow from the 
high-energy free fermion FP to the NFL FP.^^ We do not 
discuss the high-energy crossover in the present work. 

The various perturbations described by Eqs. 3 and 4 
describe very different physical processes — but the re- 
sulting crossover scale Eq. 8, has a simple form due to 
an emergent 5*0(8) symmetry of the effective NFL FP 
Hamiltonians, as discussed in the following sections. 

The main result of this paper is the NFL to FL 
crossover t matrix, given by, 

2TTiiy%a,a'a'{^,T) = Scrcr' Saa' ~ Saa,a'a'Q ( 7^ i 7^ ) (9) 

where Saa,a'a' is the scattering S matrix, which is an 
a; = and T — quantity characterizing the FL FP. For 
the 2CK model it is given by. 



S: 



2CK 

aoi.a' a' 



'(A/ 



') +»(A_B • 0„a')^aa' /A, 



(10) 



-^(7a,(T Q 



+ i(AB • Oaa')raa' I ^■ 



fll) 



The single function Q describes the crossover due to 
a generic combination of relevant perturbations in both 
2CK and 2IK models. It does not depend on details of 
the model or the particular perturbations present, except 
through the resulting crossover scale T*. Thus Q{u},T) 
is a universal function of rescaled energy tj = lu/T* and 
temperature T = T/T* . Our exact result at finite tem- 
perature is. 



\/2-k'^T 



1 _ 
2txT 






e "*■ 

dx^— — Re 

smha; 



1 1 
2^1 '2' 2'^ 



1 1 — coth X 



2t:T 



(12) 

where F is the Gamma function, and 2^1(0, &,c,z) is 
the Gauss hypergeometric function."'^ At T ~ 0, Eq. 12 
reduces^^ to the result of Ref. 26, 



g(a;,0) = -A'[-zw], 



(13) 



where K[z\ is the complete elliptic integral of the first 
kind, yielding asymptotically Q{Cj, 0) = 1 — iCj/A — 
(3w/8)2 + 0{Co^) for w < 1; and ^(^,0) = ^(tt - 
2i \og[lQCj])Cj-'^/^ + 0{Cj°) for w > 1. 

Below we consider the local density of states (spec- 
trum) fQ,o-(w,T), from which conductance can be calcu- 
lated (see Eqs. 5 and 7). It is related to the t matrix 
via Eq. 6, and is thus given exactly along the NFL to FL 
crossover by Eqs. 9-12: 



<,„(c^,r) = i-iRe 



Sa-a.aaG I ^, T 



(14) 



where the required diagonal elements of the full S matrix 
(Eqs. 10 and 11) are more simply expressed as. 



^2CK 



{-aXi + iaXg)/X = aS; 



21K 



(15) 



with (J = ±1 for spins '[ / i and a = ±1 for channel L/R 
(and we use A^ || z for simplicity). For A 



/ 



A/ 







and A^ = Xg = scattering preserves channel and spin, 
and the FL phase shift Sea then follows from Saa.aa — 
exp[2iScya]- 

These exact results for the crossover are compared with 
finite-temperature NRG calculations in Sec. VII, with ex- 
cellent agreement. 

We now examine the generic behavior of the spec- 
tral function at finite temperatures in the crossover 
regime. Although we consider explicitly L-channel spec- 
tra ta-Li^,T) in the following, note from Eq. 15 that 




Figure 3: Spectrum t^L{oJ,T) vs uj/T* for T/T* = 
10~^, 1, 10, 10^, approaching t^L = ^ from above (Ai > 0) 
or below (Ai < 0). Circles show T = result of Eq. 13. 
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Figure 4: t„L{ui = 0,T = y) and t^Liui = y^/b/a,T = 0) 
vs y/T* for Ai > 0. Common FL asymptote Eq. 16 shown 
as dot-dashed line; NFL asymptotes Eq. 17 shown as dashed 
and dotted lines. 



to-L(w,r) o tijii{uj,T) upon exchanging A^ -H- — A^ in 
the 2CK model, or Bg <-^ -Bs in the 2IK model. Also, 
i^Q(aj,T) o i^Q,(w,r) on reversing the magnetic field, 
B o — B (and in the zero-field case, a =t ^^d .^ spectra 
are of course identical). 

In Fig. 3, we take the representative case of finite chan- 
nel anisotropy A^ in the 2CK model, or finite detuning 
[K — Kc) in the 2IK model, and plot tahi'^jT) as a full 
function oiui/T* for different temperatures T/T*. Since 
only Ai acts in either case, Saa.a'a' = ^Saa'Saa' is diago- 
nal (see Eqs. 10, 11), meaning that an electron in channel 
a scatters elastically at low energies, and stays in chan- 
nel a. By Eq. (14), the spectrum t„a{LO,T) then probes 
the real part of the universal function Q because S^a.aa 
is real. 

General scaling arguments suggest that RG fiow stops 
on an energy scale given by the temperature. As seen 
from Fig. 3, this is indeed the case, with the spectrum 
tr,L{^,T) ~ tcL{0,T) essentially constant for |a;| ^ T. 
Mutatis mutandis, for T <ti T* one obtains taL{^,T) ~ 
icrL(w,0), corresponding to the T = limit consid- 
ered previously.^^ At T ~ and ui = 0, Eq. 14 yields 
ta-a{0, 0) = 5 — |ReS'o.Q,.crQ, which is determined solely by 
the S matrix and hence the phase shift associated with 
the stable FL FP. When only Ai acts, the spectrum is 
thus f 0-0(0,0) = or 1 only (with corresponding phase 
shifts or tt/2). In particular, the Kondo phase is char- 
acterized by unitarity ^0-0(0, 0) = 1, obtained in the more 
strongly coupled channel for the 2CK model, and in both 
channels for K < Kc in the 2IK model. 

In the opposite limit T ■> T* {<€-Tk) RG flow to the 
FL FP is completely cut off, and inelastic scattering^^ at 
the NFL FP results in i<,a(w,r) ~ i for all |w| < Tk- 

The generic RG picture illustrated in Fig. 1 and sup- 
ported by Fig. 3, suggests an approximate complemen- 
tarity between w and T. This is explored further in 
Fig. 4, where we compare the zero-frequency value of 



the spectrum as a function of temperature, with the 
zero-temperature spectrum as a function of frequency. 
As immediately seen, there is a strikiirg similarity. In- 
deed, a classic signature of FL physics (arising at low- 
energies/temperatures |tj|,r ^ T*) is the common 
quadratic dependence of the spectrum on both frequency 
aird temperature. 



t^a(w,r) ^ to 



«(0,0) + a( — 






(16) 



with a and b constants 0(1) that depend on the particu- 
lar model under consideration. Perturbation theory with 
respect to the FL FP in the spirit of Nozieres,"'^ yields 
a/b = 9/(77r^) for the 2IK model (see Appendix A). This 
same ratio also follows directly from the limiting behav- 
ior of the full crossover function, Eq. 12, which yields 
\a\ = 9/128 ?» 0.07 and |6| = 77r2/128 « 0.54; and as 
such provides a stringent check of our results. 

But the exact symmetry between a; and T in Eq. 16 is 
a special property of the FL FP itself, and does not in 
general apply at higher energies; although as seen from 
Fig. 4. the qualitative behavior over the full crossover is in 
fact rather similar. In the vicinity of the NFL FP (arising 
for T* <ti max(cj,T) <C Tk), Eq. 12 gives asymptotically, 



tMu^,T^O)^£'^,±{p' + Slos^}{^y-' (17a) 

1 

(17b) 



toL(a. = o,r)^^^i±/3"(^ 



T)(27-7r+log64) 

S = -^^ ~ -0.1125, 



0.5061 



with ± for Ai ^ 0; and /?' = ii±2 

(where 7 is Euler's constant), u — — 3^2 

and /3" w 0.4925. Terms of the form {ij/T*)~'^/'^ and 

(T/T*)-i/2 in Eq. 17 signal the scahng dimension 1/2 

of the relevant perturbation. Whereas such powerlaws 



occur in both the frequency and temperature depen- 
dence, additional logarithmic corrections appear in the 
frequency-dependence only. This difference can be un- 
derstood by comparing the full T = result (Eq. 13) 
with the high-T behavior captured by perturbation the- 
ory around the NFL FP (see Ref. 43 and Appendix B 2). 
The full dependence on uj and T described by Eq. 12 nat- 
urally leads to more subtle behavior when \uj\ and T are 
of the same order, as shown in Fig. 3. 

When some degree of inter-channel charge transfer is 
also present, the NFL to FL crossover is generated by 
the combination of relevant perturbations Ai and A2. 
In the 2CK model, this corresponds to finite channel 
anisotropy A^ and impurity-mediated tunneling A^, (see 
Eq. 3); while for 2IK it corresponds to finite detuning 
{K — Kc) and direct tunneling V^ji (see Eq. 4). The 
resulting behavior in the 2CK model can be simply un- 
derstood because the perturbations A^ and Aj, are re- 
lated by a 'fiavor' rotation of the bare Hamiltonian, as 
discussed further in Sec. VIA. The 2IK model does not 
possess such a flavor symmetry, although an emergent 
symmetry^^'^^ of the NFL FP Hamiltonian can be ex- 
ploited when T* <C Tk- In fact, as shown in Sec. VIC, 
this symmetry allows all the relevant perturbations in ei- 
ther 2CK or 2IK models to be simply related. ^^ implying 
the existence of a single crossover function. 

The rotation (Ai,A2) -> Ai can be used to relate 
systems where both Ai and A2 act, to those in which 
Ai alone acts (as in Fig. 3). But tcra{i^,T) probes 
the system in the original unrotated basis, and hence 
the spectra undergo a rescaling when finite A2 is in- 
cluded: t„a{uJ,T) -^ 5+1x1 [^'Tol^'^) - il- III Pai'- 

ticular, the spectral function is totally flattened the limit 
IA2I > |Ai|, with taa :i; I at both FL and NFL FPs. 
At the FL FP, electrons thus scatter elastically between 
a = L and R channels; and the corresponding S matrix 
\Saa,a'a'\ = ^crcr'(l " ^aa') IS purely Oj^^-diagoual whcu 
inter-channel charge transfer dominates (see Eqs. 10, 11). 
Thus, no crossover shows up in the spectrum or conduc- 
tance, although T* is of course finite (see Eq. 8); and the 
crossover can still appear in other physical quantities. ^^ 

When a uniform (2CK) or staggered (2IK) magnetic 
field acts (finite A^ only), Saa.aa = ii is pure imagi- 
nary (with phase shifts 5aa = i7''/4), and again we ob- 
tain <ffQ(0, 0) = i at the FL FP. However, t^rai'^, T) now 
probes the imaginary part of the universal function Q 
(see Eq. 14); and so the full spectrum along the NFL to 
FL crossover due to A^ is simply the Hilbert transform 
of the spectrum due to Ai — compare Figs. 3 and 5. 

A spectral feature in consequence appears on the in- 
termediate scale of T* for finite A|j, even though t^a = ^ 
at both NFL and FL FPs; as shown in Fig. 5. The 
existence of such a feature can be understood physi- 
cally from the impurity magnetization M ~ B^ arising 
for small applied field B^ in the 2CK model (or stag- 
gered magnetization Mg ^ B^ due to a staggered field 
in the 2IK model). Since the magnetization M{T) ex 
/^^dw /(w,T)[t^„(w,r) - t;o(w,r)] ^ is finite for 




Figure 5: t„L{uj,T) vs cu/T* for T/T* = 10"\ 1, 10, 10^ in 
the presence of finite uniform (staggered) magnetic field in the 
2CK (2IK) model, \% > 0. Spectra approach taL = \ from 
above (a =t) or below (a =^). Circles show T = Q result. 




Figure 6: Magnetization [M{T) - M{0)]/VTkT* vs T/T* 
for AI3 > and Tk/T* = 10^ 10^ 10*, lO''. NFL asymptote 
Eq. 20b shown as dashed line; FL asymptote Eq. 20a shown 
as dotted line in the inset. 



finite applied field, t^a{^,T) 7^ i^a(a;,r). A 'pocket' 
thus opens between the 'up' and 'down' spin spectra at 
\uj\ ^ T* , whose area is proportional to the magnetiza- 
tion. In particular, the temperature-dependence of the 
magnetization can be extracted from the universal func- 
tion, Eq. 12, viz 



M{T) oc 



Tk 



duj f{oj,T) Img ( — ,— 



(18) 



valid for small perturbations A|j, such that T* ^ Tk 
as usual (the high-frequency cutoff |w| ^ Tk then being 
justified since i^a(w,r) ~ tia{^,T) for \uj\ > Tk, as 
confirmed directly from NRG). 

Using Eq. 13, the zero-temperature magnetization de- 
pends on Tk via. 



M(0) = c^T^ \og (^^) 



(19) 



and since T* - {B'^f/TK (see Eq. 8), it follows that 
M(0) - Bnog{TK/B''). The fuh temperature depen- 
dence M{T) vs T is shown in Fig. 6, demonstrating scal- 
ing collapse for different Kondo scales Tk- The asymp- 
totic behavior in the vicinity of the FL and NFL FPs 
follows as, 



VWT 



^ 'M{T) - M(0)] ^/^«^ 






T/T*»l ( T 

-clog!—), 



X™p(r) ~ \og{TK/T), 



(20a) 

(20b) 

yielding in particular M{T) - B^ log(TK/T) when T* < 
T ^ Tjf . This asymptotic behavior can again be under- 
stood from perturbation theory around the FL and NFL 
FPs. Furthermore, since X\mp{T) = limB._,o Af (^)/B^ 
when T* <^ B^ <s:T <t.TK, one obtains 

(21) 

for the uniform (staggered) magnetic susceptibility of the 
2CK (2IK) model. This diverging susceptibility is a clas- 
sic signature of the NFL FP, known for example from 
the Bethe ansatz solution of the 2CK model,"''* or from 
CFT for the 2IK model. ^'^ However, at lower tempera- 
tures T < B^, M{T)/B'' does not correspond to the 
magnetic susceptibility: here the NFL to FL crossover 
itself is being probed. Indeed, for T ^ T* < B^, we 
obtain a quadratic (T/T*)'^ temperature-dependence of 
magnetization, Eq. 20a, characteristic of the FL FP. 

Finally, we turn to conductance G"{V,T), obtained 
from the spectrum tcra{'-^,T) by combining Eqs. 7, 12 
and 14. It follows as 



1 






'^ ^ 2 (87rr)3/2r(l 



1^ ■ 

2ttT' 



z^ 



Im< Saa., 



dx I{V,T, x)Rc 



2F1 



1 1 — coth X 



2ttT 



(22) 

where the integral over lo can be evaluated using contour 
methods. 



IiV,T,x) 



du) 



exp 



/ ixu \ 



sech 



2T 



sinh(a::) tanh (^) 



--2Trif csch^(a;) sech^ (^) 



cosh(x) 



(23) 



^-P(#)(l + Tsinli(f 



with rescaled V = V/T* , T = T/T*, u = uo/T* as before. 
In particular at zero-bias. 



/(y = o,r,x) = ^^ 



cosh^ (f ) 



(24) 



G^(V,T) 



LO 




Figure 7: Conductance G^{V,T) = G^{V,T)/{2e'^h-''G^) vs 
bias V and temperature T for Ai > 0. Black lines connect 
regions of constant conductance. Light colors correspond to 
high conductance near the FL FP; dark colors correspond to 
lower conductance near the NFL FP. 



A color plot of conductance G"{V,T) along the NFL 
to FL crossover is shown in Fig. 7, for the representative 
case of Ai > 0. The black lines connect regions of equal 
conductance. In the FL regime ]/, T ^ T* , the quadratic 
form of Eq. 16 thus yields a simple ellipse; while in the 
NFL regime there is a pronounced V-T anisotropy. The 
detailed behavior of Eq. 22 is seen on taking cuts through 
Fig. 7 at constant V and T, as shown in Fig. 8. 

Returning to the 2CK quantum dot system of Ref. 24, 
we comment now on the possible strength of symmetry- 
breaking perturbations present in the experiment. Due to 
the Coulomb blockade physics of the quantum box, inter- 
channel charge-transfer was effectively suppressed. In the 
absence of a magnetic field, the dominant perturbation is 
thus channel anisotropy, A^ (see Eq. 3). The experiment 
showed^^ 2CK scaling of conductance around Tk , but no 
FL crossover at T*. This implies T* ^ T <^ Tk — see 
for example Fig. 3 for T/T* = 100, which shows little sign 
of the crossover. Since T /Tk ~ 0.1 in the experiment and 
taking T/T* > 100, from Eq. 8 and Table I it follows that 
cii^Az could be at most ~ 0.03 (ci = C(l) depends on 
details of the model/device setup). The observed 2CK 
physics is thus an impressive testament to the tunability 
and control available in such quantum dot devices. 

Having presented our main results and discussed their 
physical implications, we turn in the following sections 
to the formal derivation. 
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Figure 8: Conductance G^(V,r) = G^{y,T)l{2e^h-^Go) 
arising for finite Ai . Upper panel: vs temperature T/T* for 
V/T* = 10"\l,10, 10^ approaching GJ = f from above 
(Ai > 0) or below (Ai < 0). Circles show the exact zero-bias 
result. Lower panel: vs bias V/T* for T/T* = 10~\ 1, 10, 10^ 
Diamonds show zero-temperature result. 



tron operators can be Fourier transformed and expanded 
near the Fermi points, focusing on states within width 
2D <C vpkp around ei?, 



*.a(x) - ^e'^'^-V'fc 



(25) 



kp+D/vE 
k—kf—D/vp 
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We thus define left and right movers, 

D/VF 

k=-D/vF 



(26) 



defined for x > 0, with x the distance from the im- 
purities located at the 'boundary' a; = 0. With a 
boundary condition 'ipr,aa{Q) = ''Pi,aa{'^) we introduce 
a single left-moving chiral Dirac fermion for all x, 

i^Ga{.x) = 9{x)ll}l^cra{x) + 6'(-x)V'r-,aa ("X). SinCC the 

new operators satisfy the usual fermionic anticommuta- 

tion relations {^la{xi),i^„,a'{x2)} = Saa'Saa'S{xi -X2), 

the free Hamiltonian Hq can be expressed as a chiral 
Hamiltonian."*^ 



H^^vp^l dx^l^id^i:,^. (27) 



III. EXACT T = Q CROSSOVER GREEN 
FUNCTION IN THE 2IK MODEL 

Our goal is an exact expression for the NFL to FL 
crossover t matrix, which is related to the electron Green 
function. In Ref. 26, we calculated the crossover at T = 
0; further details of that calculation are presented here, 
providing as they do the necessary foundations for our 
generalization of the results to finite temperature. 



A. Fixed point Hamiltonians and Green functions 

The structure of the NFL fixed point Hamiltonian of 
the 2IK model allows for an elegant description of the 
NFL to FL crossover. ^^ Before presenting that derivation, 
we discuss first some relevant preliminaries which will be 
of later use. In particular, we consider now the represen- 
tation of the fixed point Hamiltonians within CFT, and 
the structure of the corresponding Green functions. 

Our starting point is a description of the free conduc- 
tion electron Hamiltonian Hq in terms of chiral Dirac 
fermions. A ID quadratic dispersion relation e(fc) = 
Efc = k^ /2m— ep can be linearized near the Fermi points 
k = ±kp, €k — ±vp{k ^ kp) (with vp the Fermi ve- 
locity). This is the standard case^ and applies to ar- 
bitrary dimension within the assumption that the bare 
density of states is fiat at low energies.'*'^ Conduction elec- 



Hereafter, we set vp = 1 [and the Fermi level density of 
states is then 1/ = l/{2Trvp) = l/(27r)]. 

The free electron Green function then follows as,^^ 



(^.a(r,Xi)^L'(0,^2))o--T^ 






&i-n[^{T + ixi - 1x2)] (28) 

5aa'5aa'G'^{T,Xi - .T2), 



where r is imaginary time. The corresponding Matsub- 
ara Green function is then defined by the transformations 



G 



V,xi-x.) = i^e— ° 



/3 



Q {Xl ~ X2,iUJn) 



P/2 



(29) 



g"ix,-X2,itUn)^ I rfTe'""^G°(r,xi-X2), 

-13/2 



with Matsubara frequencies aj„ = ^(1 + 2n) for integer 
n, and where /? = l/fc^T is inverse temperature. Direct 
evaluation of Eq. 29 using Eq. 28 then yields simply, 



Sacr'Saa'[0iuJn)dix2 ~ X\) - 9{-UJn)0{xi - X2)]. 



(30) 



Of course, the interesting behavior arises when cou- 
pling to the impurities is switched on. As usual, ^ the full 
Green function is related to the t matrix via 



g{xi,X2,iuJn) = Q°{xi -X2,iUJn) 
+ g°{xi,iuJn)T{iL0n)G°i-X2,iuJn), 



(31) 



where Q^ Q and T are 4x4 matrices with indices taking 
the values aa =t L, 4- i, t -Ri i ^- In particular, it should 
be noted that the t matrix is local. Eqs. 30 and 31 also 
imply that if .ti and xi have equal sign then the full Green 
function reduces to the free Green function, reflecting 
the chiral nature of the Dirac fermion, Eq. 27. Since 
all such fermions are now left-moving, information about 
scattering from the impurities located at the boundary 
X = — and hence the t matrix — is obtained from the 
full Green function with X\ and X2 located on opposite 
sides of the boundary. 

The free Green function Eq. 28 is naturally obtained at 
FPs where the free boundary condition pertains. But at 
any conformally invariant FP, the powerful machinery of 
boundary GET gives nonperturbative information about 
the Green function. Specifically, when the boundary con- 
dition is obtained from fusing with some primary field a, 
then correlation functions are given generically by"^^ 



{O^{T,tX,)Oli0,tX2)) 






\t -f ixi — ix2 



2d' 



(32) 



where d is the scaling dimension of the primary field O, 
and 5*" are elements of the modular S matrix. '^■^ Using 
the conformal mapping from the plane to the cylinder 
with circumference /?, one obtains^^ the generalization 
to finite temperature T = j3~^, 



{0^{T,lX^)0\{Q,lX2)) 






(^sin[f(r + *xi-zx2)])2'i 

(33) 
In the present context, we are interested in the electron 
Green function at the conformally invariant free fermion, 
NFL and EL EPs of the 2IK model. With O = ifj^a the 
d = 1/2 fermion field, the full EP Green functions takes 
the form,^^ 



Gaa,a>a>{T,Xi - X2) 



1 TT C 

2ir p'^<ya,a'a' 



sin[|(T + ixi - 1x2)]' 



(34) 



where S can be understood as the one-particle to one- 
particle scattering matrix, and can be calculated from 
the modular S matrix in the case of boundary conditions 
obtained by fusion. '^^ Thus, the effective FP theory is 
identical to that of the free theory discussed above, but 
with a modified boundary condition that determines the 
scattering matrix 5'. 

Ghoosing xi > and X2 < in Eq. 31 yields i7~ = 
1 — 5^, or equivalently ta-a = 5 (1 ~ ReS'o-a.cra), where the 
above convention v = l/(27r) was used. At a EL FP, the 
scattering matrix is unitary,^ 5'^ 5' — 1, and as such de- 
scribes purely elastic scattering. By contrast, at the NFL 
FP of the 2IK model, it has been shown^^ that S = 0, 
which implies fully inelastic scattering: a single electron 
sent in to scatter off the impurities decays completely 
into collective excitations, and no single-particle state 
emerges. Such behavior is manifest by a half-unitarity 
spectrum, tcra = k- 



However, along a crossover between FPs, the Green 
function does not in general take the form of Eq. 34. 



B. S0(8) Majorana fermion representation 

Further insight into the FPs of the 2IK model is pro- 
vided by a representation in terms of Majorana fermions 
(MEs). Considering again the free theory described by 
Hq, four nonlocal fermions can be defined by Abelian 
bosonization and refermionization^^'^'''^^ of the four orig- 
inal Dirac fermions ip^a with spin a =t: i ^-^id channel 
index a ^ L, R. 8 MEs are then obtained by taking the 
real and imaginary part of each. 



Specifically, four bosonic fields (j)^a are defined, viz 



V'a 



-^ (70.^ 



(35) 



where F^ra are Klein factors.^® Linear combinations of 
these bosonic fields are then used to construct new fields, 

{^e, </>.>/, 0X} = i Y. -^-{l' (-1)"^'' (-l)"+\ (-1)"+"}, 
act 

(36) 
where a =^, 4,= 1, 2 and a= L, R= 1, 2. Refermionizing 
yields. 



IpA ^ FaC 



-itpA 



(A 



J,X), 



(37) 



where the Klein factors Fa are related to Fg-a as de- 
scribed in Appendix C. Thus, four new fermionic species 
ipA are defined, with A = c,s,f,X corresponding to 
'charge', 'spin', 'flavor' and 'spin-flavor'. The real and 
imaginary parts of each, 

A _ i^A + i'A A _ V'a-V'A ,„„^ 

Xl - 7= , X2 — 7^. ' l-J^J 



V2 



V2i 



fulfill the Majorana property (x^*)^ = xf and sat- 
isfy separately the fermionic anticommutation relation 
{{xt)H^)^xf{^')} ~ SijSab^X'-x'), and so are referred 
to as MEs. 

The free fermion FP Hamiltonian (corresponding to 
Eqs. 1 and 2 with J = and 6H = 0) then follows as 



Hfp = ^ 



dx x{x) ■ d^xix), 



(39) 



where x= {X2 ^x{,xLxl,X2^Xi ,Xi,X2}^ and the scat- 
tering states are defined by the trivial boundary condi- 
tion xi^) = x{~^)\x^o- The fixed point thus possesses 
a large S0{8) symmetry in terms of these MEs. 

Likewise, the EL FP Hamiltonian (in which the im- 
purity degrees of freedom are quenched) is similarly de- 
scribed by Eq. 39. The corresponding boundary condi- 
tion is encoded in the single-particle EL scattering S ma- 
trix Sa-a.tr'a'- Although it depends on the specific pertur- 
bations generating the crossover, the boundary condition 
is thus trivial at the EL FP. In particular, finite detuning 
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K > Kc results in x(a;) = x{^^')\x^o, as obtained at the full Green function coming from each of the sectors: 
free fermion FP. 

7 

The remarkable fact is that the NFL FP Hamiltonian r^NFL 

also takes the form of Eq. 39; and its nontrivial boundary 
condition^"^ is again simple in terms of the MFs. It can be 
accounted for by defining the scattering states x^(x) = 
-X2i-x)\x^o, and xf{x) = xf{-x)\x^o for (j, A) ^ 
(2,X). Thus, 7 of the 8 MFs are described by the free 
theory at the NFL FP. 



C. Bose-Ising decomposition 

The FP Hamiltonian Eq. 39 describes a higher SO{8) 
symmetry than is present in the original Hamiltonian. 
The explicit symmetries of the 2IK model also allow a 
separation of the theory into different symmetry sectors. 
Specifically, the 5't/(2)ixS't/(2)ixS'[/(2)2xZ2 symmetry 
sectors comprise a Bose-Ising representation,^'^ describ- 
ing a coset construction of three Wess-Zumino-Witten 
(WZW) nonlinear a models, together with a Z2 Ising 
model. The two SU{2)i theories with central charge 
c = 1 correspond to conserved charge in the left and 
right channels. The SU{2)2 theory with c ~ 3/2 corre- 
sponds to conserved total spin. Finally, the Ising model 
^2 is a c = 1/2 theory corresponding to a single MF. This 
non-Abelian representation is connected with the 8 MFs, 
as discussed in Ref. 47. The symmetry 'currents' of those 
sectors, such as the spin current J{x) — ip''^ {x)^ip{x), are 
represented quadratically in terms of MFs as described in 
Appendix C. Specifically the SU{2)i x SU{2)i charge 
currents in left and right channels are represented in 
terms of 4 MFs {x{, xi) and {xl-, X2}; while the 51/(2)2 
spin current J is represented in terms of three MFs, 
Xs = {Xi,X2;Xi'^}- The Z2 theory corresponds to the 
remaining MF, X2^- 

The important implication for our purposes is that 
the Green function can be factorized into pieces coming 
from different sectors associated with the various MFs. 
We now exploit the above Bose-Ising construction,^'^ in 
terms of which the fermion field can be expressed by the 
bosonization formula, 



1paa[x) CX [ha\l{x)g„[x)(7L{x). 



(40) 



Here, the dimension d = 1/2 fermion field has been de- 
composed into a dimension dh = 1/4 factor [ha\i repre- 
senting the a ^ L,R spin-^ primary field of the SU{2)i 
charge theories, a dimension dg ~ 3/16 factor {g^) repre- 
senting the spin-2 primary field of the S'L'^(2)2 spin the- 
ory, and the dimension da = 1/16 factor ai^ originating 
from the Ising sector. The subscript L emphasizes that 
(Tl is only the left-moving chiral component of the full 
spin operator a of the Ising sector, arising here because 
■00-0! is the chiral left-moving fermion field. Since the NFL 
FP of the 2IK model is conformally invariant, ^^ we may 
use Eqs. 40 and 33 to determine the contribution to the 



G'^Q^i''a'(T, Xl - X2) (X baa'&aa' [G°(t,. Ti - X2^ ^ I/^y\ 

X ((Ti,(t, ia;i)cri (0,1x2)) ■ 

The free boundary condition in the charge and spin sec- 
tors yields the first factor, corresponding to the free 
Green function Eq. 28 but with power 2(d^ -I- dg) = | 
arising because 7 of the 8 MFs are associated with 
these sectors. The NFL boundary condition is expressed 
in terms of fusing with the da ~ 1/16 field ul from 
the Ising sector in Eq. 40. The second factor thus 
comes from the remaining Ising sector. At the NFL 
FP it follows from Eq. 33 that (crL(r, zxi)crL(0, 1x2)) = 



1/16/q1/16\ j_ 

T^l [G0(t,Xi-X2)]8. 






Since the modular S 



matrix for fusion with the a ~ 1/16 Ising operator has a 
vanishing element S-^L^ = 0, the entire Green function 
thus vanishes at the NFL FP;^'^ consistent with Eq. 34 

with Saa.a'a' = 0. 

In summary, the nontrivial boundary condition at the 
NFL FP affects only the Ising sector of the 2IK model. 
The function ((TL(r,ixi)crL(0,zx2)) in Eq. 41 is a quan- 
tity pertaining to the Z2 Ising model and contains the 
nontrivial physics; while the spin and charge sectors are 
simply spectators. In the next section we exploit Eq. 41 
and a connection between the 2IK model and a classical 
Ising model^^ to determine the Green function along the 
crossover from the NFL FP to the FL FP. 



D. Crossover due to K j^ Kc in the 2IK model 

As highlighted above, the NFL FP of the 2IK model 
is described by the free theory in all sectors except the 
Ising sector, which takes a modified boundary condition. 
When the small perturbation Ai is included (correspond- 
ing to detuning K 7^ Kc), the NFL FP is destabilized. 
The effective Hamiltonian in the vicinity of the NFL FP 
is then Hqcp = Hfp[x\ + SHqcp, with Hpp the FP 
Hamiltonian itself, given in Eq. 39 and parametrized in 
terms of the scattering states x which encode the bound- 
ary condition. The correction is given by. 



SHqcp oc iAiX2 (0)a 



(42) 



where a is a local MF involving impurity spin operators^^ 
which satisfies a^ ~ 1 and anticommutes with all other 

MFs, x^ 

The perturbation Ai thus acts only in the Ising sector. 
Indeed, the difference in boundary conditions between 
the NFL and FL FPs is also confined to the Ising sector 
(all other sectors have free boundary conditions at both 
FPs). The entire crossover from NFL to FL FP thus 
occurs completely within the Ising sector because the Ai 
perturbation does not spoil the decoupling of the sectors, 
which becomes exact^^ in the limit T* <C Tk- The other 
sectors then act as spectators along the crossover. 
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This implies a generalization of Eq. 41 to the full 
crossover. One can still interpret the first factor in Eq. 41 
with power 7/8 as the product of autocorrelation func- 
tions of 7 spin fields which undergo no change in the 
boundary condition. But the a field originating from the 
Z2 sector in the SU{2)i x SU{2)i x SU{2)2 x Z2 con- 
struction has flowing boundary condition (which is not 
conformally invariant) . 

We consider first the equal-time Green function with 
xi = —xi positioned symmetrically on either side of the 
boundary: 

7 (43) 

oc [G°(0,2.t)]8 X {GL{T,ix)aL{T,~ix)), 

where the factor (ctl(t, ia;)CTL(T, —ix)) now describes the 
crossover in the Ising sector in terms of the chiral (left- 
moving) Ising magnetization operator cri(r, ice). 

We now use standard boundary CFT methods'*^ to re- 
late the two-point chiral function (cri(r, ix)crL(r, — ix)) 
living in the full plane [see Fig. 9 (b)] to the product 
of chiral holomorphic (7]^(t + ix) and antiholomorphic 
(Tfl(r — ix) operators living in the halfplane with a bound- 
ary. But bulk operators in CFT may be expressed as 
cr(T, a;) = (Jl{t + ix)a^{T — ix), and so the desired corre- 
lator is simply (crL{T, ix)(JL{T, —ix)) = {(j(t,x)) in terms 
of the bulk Ising magnetization operator, evaluated at a 
distance x from the boundary [see Fig. 9 (a)[. Finally, 
we note that cr(T, x) = <y{x) is independent of t due to 
translational invariance along the boundary. The Green 
function along the NFL to FL crossover then follows as, 

1 

G<,„,,-„'(0,x,-x) ^5,,,5,,^, (^y [G°(0,2x)]8 (a(x-)), 



(a) 



-^(TOi^a' ol' y^) 1 



(44) 



(with the factor (Svri) ^/^ required for normalization). 
At long distances where the Green function describes the 
FL FP, Eq. 44 implies {a(x)) = (2/x)i/« as T -^ 0. 

We now exploit a connection'^'^ between the 2IK model 
at criticality with a simpler Ising model to obtain {(y(x)) 
and hence the exact Green function Gaa,a'a' {x) along the 
NFL to FL crossover. 

When a small magnetic field h is applied to the bound- 
ary of a classical Ising model on a semi-infinite plane, 
the local magnetization shows a crossover''^''^^ as a func- 
tion of distance from the boundary. The crossover in this 
boundary Ising model (BIM) can be understood as an RG 
flow'^^''^^ from an unstable FP at short distances (with 
free boundary condition h = {)), to a stable FP at large 
distances (with fixed boundary condition h = ±00). The 
universal crossover is characterized by an energy scale 
T* = Airh'^ (or a corresponding lengthscale £,* (x 1/T*). 

Importantly, it was shown in Ref. 23 that the RG flow 
in the BIM due to small h at the critical temperature is 
identical to the RG flow in the 2IK model due to small 



X 

i 



(b) 



• Z=IX 



^z*=-ix 



(c) 



Z=X+IX ^ 



Z = Tpx^ 



Figure 9: A one-point function of a bulk operator evaluated 
at distance x from the boundary in (a) is mapped to the two 
point function of the associated chiral fields in the absence 
of a boundary sliown in (b) and evaluated at image positions 
with respect to the line a; = 0. (c) Generalization of the two 
point function away from image points. 



detuning A' 7^ iCc at T = 0. As such, the NFL and 
FL FPs of the 2IK model can be understood in terms of 
the BIM FPs with free and fixed boundary conditions. 
The crossover energy scale in the 2IK model can then be 
identified as 



T* ^ ^^ — = Xl 



Tk 



Airh^ 



2IK 

BIM 



(45) 



with ci = 0(1) as in Table I. 

In Ref. 33 Chatterjee and Zamolodchikov derived an 
exact expression for the Ising magnetization ((t(x)) on 
the semi-infinite plane geometry in the continuum limit. 
Their result^^ is 

{^{x))cz = ±{2/x)^/'^V8h^ e'^^''''=Koi'iTrh^x), (46) 

with Kq the modified Bessel function of the second kind; 
and ± for ft, ^ 0. We take now h > (corresponding 
to K > Kc) for concreteness. Note that Eq. 46 yields 
asymptotically (o'(x)) = {2/x)^/^ at long distances x — >■ 
00, consistent with the normalization of Eq. 44. 

We now show that the analyticity of the Green func- 
tion and the local nature of the Kondo interaction im- 
plies a generalization of Gaa,(y'a'{x) in terms of spa- 
tial coordinate x [Fig. 9 (b)], to Grja,<T'a'{zi — 2^2) = 
('i/'cra(-2i)^CT'Q'(-Z2)) in terms of general complex coordi- 
nates zi ^ T -\~ ixi and Z2 = 1x2 [Fig. 9 (c)]. Using the 
free chiral Green function Eq. 30 in the definition of the 
t matrix Eq. 31, one obtains 

Gaa,a'a'{xi,X2,iuJn) = -S^a' Saa'^i-i^n)^ 

e^-^^'-^'^[i + %o..a'a'{iOJn)], 

(47) 

where xi > and X2 < as before, and the t matrix is 
local in space. 
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The Matsubara transform then yields 



^(Ta,{T'a' \^5 ^1 T -^2 ) 



E^(-^») 



X 






— ^(Ta.cr'a' 



2i 



(48) 



Thus, the Green function only depends on (zi — 22)- 
This is a somewhat counter-intuitive resuh, because the 
boundary breaks the translational invariance along the 
spatial coordinate x. In Appendix B 1 we give an alter- 
native proof, showing that Eq. 48 holds to all orders in 
perturbation theory around the NFL FP. 

Comparing Eq. 48 and Eq. 44, it follows that 

Gaa,a'a'iO,X,-x) = Gcya,a' a' (t, Xi, X2) whcu X ^ (zi ~ 

Z2)/{2i). Employing this substitution in Eqs. 44 and 46, 
then taking the Matsubara transform, we obtain 



ycra,(T'a' V*^l ; -^2 ; ^^n J — ^aa'^ac 



dT 



VSh^x 



i-Kh X 



Airix 

2 



Koiinh^x), 



(49) 



where we have used /3 — > oo as appropriate for Eq. 46, 
and so a;„ ~ ^(1 + 2n) is continuous. Setting xi = 0+, 
X2 = 0^, we define the infinitesimal S = xi ~ X2 > 0, 
such that X — ^-^^- With the integral representation 

of the Bessel function Kniz') ~ e"^ L dk f (for 

^ ' •'0 y/k{k+l) ^ 

Re z > 0), we obtain 

g{ujn + 4:TTh'^k) 



Try J -'0 



(50) 



with the integral over r evaluated by contour methods 

(>oo-\-iS 



3(^) 



dT^ = e{-z)J—. (51) 

oo+iS V^ V Z 



For negative Matsubara frequencies w„ < 0, the Green 
function then follows as 



Gaa,cr'a'{0^,0 , iw„) 



dk 







Vkik + l)i~i^~k) 



■Ki V 47rft,2 / 



(52) 



where the k integral has been expressed more simply 
in the last line in terms of the complete elliptic in- 

r7r/2 



tegral of the first kind, K{z) = L —i= 
ri dt 

•^0 y/{l-t^)(l-zt^)' 



Since K{z) has a branch cut discontinuity in the com- 
plex z-plane running from 1 to c», the analytic contin- 
uation to real frequencies, iw„ -^ cu + iO+ can be per- 
formed without crossing any singularity. Thus, if one 
has Ga-a,a'a'ix) as an analytic function of spatial coordi- 
nate X (as in Eq. 44), then a full knowledge of both space 
and time dependences of the Green function is implied 
by analytic continuation. 

Using Eq. 31, one recovers our earlier result^^ for the 
T = crossover t matrix in the 2IK model due to per- 
turbation K y^ Kc, 



2TTiiy'Ta 



(TQ.cr'Q 



(w,r = 0) = 5,,.(5„„, lT-K\ 



-iLu/T* 



(53) 
in terms of T* = Anh'^ (and with v = l/(27r) as before). 
Here ^ is used for K ^ Kc, corresponding to local singlet 
or Kondo screened phases of the 2IK model, respectively. 
Eq. 53 is thus equivalent to Eqs. 9 and 13 with the scat- 
tering S matrix Sg-a.a'a' = ^^a-a'^aa'- 

In the next section, we generalize these results to finite 
temperature. 



IV. DERIVATION AT FINITE 
TEMPERATURES 



Our starting point for the derivation of the finite- 
temperature crossover Green function is Eq. 44. 
Gaa.ij'a'{x) thus follows from the Ising magnetization 
{a{x,l3)) evaluated at temperature T = (3^^ . 

In Ref. 43 we considered the magnetization at tempera- 
ture l3~^ and distance x from the boundary in a quantum 
ID transverse field Ising critical chain, with a magnetic 
field h applied to the first spin at the point boundary. It 
is given by 



{a{x,P))=f{Ph^)x{aix,P)) 



LLS, 



(54) 



with 



^(^"^' = ^^ wi^y '''' 



and where 



(a(x,/3)) 



LLS 



^sinh2^ 



/l 1 i_coth2^\ 
2F,[^-,^jA + 2ph^ T^j' 

(56) 

is the result of Leclair, Lesage and Saleur in Ref. 34, 
who generalized the T = result of Chatterjee and 
Zamolodchikov'^'^ for the semi-infinite plane geometry 
(Eq. 46) to the geometry of a semi-infinite cylinder with 
circumference /3, and magnetic field h > applied to the 



13 



circular boundary. The latter is equivalent to the quan- 
tum Ising chain with transverse field. We showed*'^ that 
while {a{x,f3))LLS gives the full and highly nontrivial x 
dependence of {a{x, P)), it misses the multiplicative scal- 
ing function of the variable /3/i^, given in Eq. 55. 

In the low-temperature limit /3 — > oo, one obtains 
fiPh'^) -^ 1 and 

(l 1 l-coth^\ ,,^^ 

{Sh^xf'^e^^'^'^K^i'iTTh^x), 

such that Eqs. 54, 56 reduce as they must to Eq. 46. In 
the limit /i — > oo, one recovers the fixed boundary con- 
dition, describing the FL fixed point. Again f{fih^) -^ 1 
but Eq. 56 reduces now to 



{a{x,P)) 



LLS 



^sinh^^ 



(58) 



A. Ambiguities in analytic continuation at finite-T 

The quantity of interest is of course the t matrix 
TiiiOn), related via Eq. 31 to the Matsubara Green func- 
tion, itself obtained from Eq. 60 via 



G-tL,tL{xi,X2,iu}n) = I dr 6*'^"'^ G- 



13/2 
-13/2 



t-L.ti 



^1 - Z2 

2i 



(61) 

with zi = T + ixi and Z2 = ix2 as usual. Using the 
integral representation of the hypergeometric function, 

(62) 
we then obtain 



ytL,tL{0 ,0 ,tUJn) = 3/2^ / at y= An, 

(63) 

where 



This limiting behavior can also be obtained using the 
conformal mapping from the semi-infinite plane geome- 
try with boundary Imz = 0, to the semi-infinite cylin- 
der geometry parametrized by Rez' e (— /3/2,/3/2), with 
boundary Imz' = 0, 



tan 



TTZ 

T 



(59) 



On the semi-infinite plane, the Ising magnetization in the 
limit /i — >■ oo is known^^ to decay as {a^x)) = {2/xy^^, 
yielding precisely Eq. 58 on the semi-infinite cylinder. 

Combining Eqs. 44, 54-56 we obtain the crossover 
Green function at finite temperature. 



G. 



(TOL.a' a. 



{J^j — ^aa' ^cxa' 



Vrn^ 



r[i + 2/3/i2 



1 1 



r[l + 2/3/i2] 2^1 I 2' 2'^ 



2^isinh[^] 

1 - coth 2f^ ^ 
2ph\ ^ 



(60) 



For h — ^ oo, this gives correctly Gcra,a'a'{0,x, —x) — 
Gf^^,„,(0,x,-a;) = l/(2/3isinh[2|^]), as expected from 
the boundary CFT result for the FL FP Green function, 

Eq. 34 (with S^a,aa = 1)- 

Of course, the function f{ph?) does become impor- 
tant when considering the behavior of the Green func- 
tion over the entire range of /J/i^ (or equivalently, T/T*). 
In particular, at high temperatures /3 — J> (and finite 
h), one obtains f{/3h'^) -> 0. Thus Gaa.a'a'{0,x,—x) = 
G^^^,^,{0,x,~x) — 0, again correctly recovering the 
known boundary CFT result^^ for the NFL FP Green 
function, Eq. 34 (with Saa,cr'a' = 0). The factor /(/3/i^) is 
indeed necessary to cancel the unphysical logarithmic di- 
vergence of {(j{x, P)) LLS as /3 — )■ (see also Appendix A). 



An = dr e fi 



(2n+l) s\nh{f{S-tT)) 



(64) 



for xi ~ 0"*" and X2 = 0~ such that S — xi — X2 > 0; 
and with n — '^^ ~ \ ^ negative integer. Using contour 
integration, it can be shown that 



^mii)n+iJ}_Jl 



-ri-1 



An = 8^^/^(-l) 



T{~n)T{l + n) 



,1 3 1 

2^1 -,1 + n, - + n. 



2 ' l-< 



(65) 



However, the naive substitution 



uj (or n — >■ 



^ — ^) is problematic here and leads to unphysical 
divergences. Indeed, such analytic continuation always 
involves ambiguities due to the fact that (— 1)^" = 1 
on the integers, but it becomes —e^^ upon analytic 
continuation. Thus, it is hard to find the function 
Qi;L,^l{Q~^ 1 0~, iLOn) which gives the physical analytic con- 
tinuation. 



B. Finite-T Green function from 
Friedel oscillations 



Eq. 60 describes the chiral electron Green function 
Gaa,a'a'{x) = {^Jaaix)^]^, ^,{—x)) . The information con- 
tained in such Green functions is directly linked to the 
physical density oscillations around impurity (Friedel os- 
cillations), which in turn can be calculated from the t 
matrix, 7~(a;).'^'^'^^'^° Indeed, in Ref. 35 the real-space 
densities and hence Green function Gaa,iya{x) for the 
NFL to FL crossover in the 2CK model was explicitly 
calculated at T = using the exact t matrix announced 
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in Ref. 26. It was also highlighted^^ that far from the im- 
purity, the integral transformation relating the t matrix 
to the Friedel oscillations can be inverted. 

In this section we exploit these connections to calculate 
T{uj) directly from the density oscillations described by 
Gaa,a-aix), and thus circumvent the need for problematic 
analytic continuation. 

For simplicity we restrict ourselves to ID in this sec- 
tion; although we note that the resulting t matrix is 
general, because at low energies the standard flat band 
situation of most interest is recovered. The density of 
the ID fermion field \E'(x) at position x is given by 
p{x) = {'i^x)'^ {x)) . Expanding around the left {I) and 
right (r) Fermi points at low energies using Eq. 26, one 
obtains 



*(.T)=^,(x)e*'=^"+Vi(a;)e 



-ik^x 



(66) 



with the oscillating part of the density following as 



Posc{x) = ie2"=^"(V'/(x)V'r(x)-^.(a:)^/(a;))+H.c. 



(67) 



In the presence of particle-hole symmetry, p{x) = 1/2 for 
all X (with lattice spacing set to unity) ; and there are no 
density oscillations, Posc{x) = 0. However, introduction 
of potential scattering breaks particle- hole symmetry and 
leads generically to real-space density oscillations, which 
contain information about the t matrix. Such potential 
scattering produces a phase-shift 5p at the Fermi energy, 
independent of the underlying Kondo physics, but which 
does modify the boundary condition at x = 0, accord- 
ing to '0^(0) = e^*''^'0i(O). As before, we define a chiral 
left-moving field on the infinite line, but now take into 
account this change in the boundary condition: 



^{x) = Mx)0{x) + e-^'''"'iPr{-x)e{-x) 



(68) 



Using the imaginary time ordering of the chiral Green 
function defined in Eq. 44, we have G{~x + iO+) = 
{tp{—x)i}/^{x)) and G{—x — iQ^) — —{tp'^ {x)^{—x)), from 
which it follows that 



Pose (a;) 



^ 2ikFX-\-2iSp 
2^ 



+ Gaa.aa{-X + iO+) 



H.C.. 
(69) 



The oscillating part of the density given by Eq. 69 can 
also be obtained'^'''''*^''''° from the t matrix. Generalizing 
to finite temperatures, we have 



Posc(X) 



duj fiLu,T)x 
Im{[G°{u,x)]^%aM^,^p)} , 



(70) 



where f{Lu,T) is the Fermi function, G'^{lu,x) is the free 
Green function Eq. 28 as a function of real frequency w; 
and 7^Q,o-a(w, (5p) is the scattering t matrix, defined in 



the presence of the potential scattering. As shown in 
Ref. 50, at low energies 



2TTvTac 



,{uj,Sf 



2iSp 



2lTV%c 



.(w) 



(71) 



in terms of the desired t matrix defined without the po- 
tential scattering. Indeed, far from the impurity one ob- 
tains asymptotically'''' 

[G°(W, X)f = ^J^^2^k^x+2^ux/v^ _ (73) 

The oscillating part of the density can then be ex- 
pressed as 



Pose \X) 



Attv 



^— f doj[l-2f{uj,T)] 



_2ikFX-\-2i5p--\-2iijjxlvf 



i%i 



1 

2inj 



H.c. 



(73) 



Comparing now to Eq. 69 and inverting the Fourier 
transform by operating on the resulting equations with 
r°° dx^-2%uix ^^ obtain 

J — oo 7r ' 



2'KivTaa,cra{^) =1 



Anvvp 



dx e^^"^/"^ X 



tanh(^) J-oo 

[Gcra,aa{x - iO + ) -|- Gaa,cra{x + iO+)], 

(74) 

where Gaa,<Ta{x) is given as an analytic function in 
Eq. 60. We now set Ui? = 1 and 27r!^ = 1 as before. 
Note that the hypergeometric function 2F1 (a, b, c, z) has 
a branch cut discontinuity in the complex z plane run- 
ning from 1 to 00. The discontinuity occurs only in the 
imaginary part of the function, with Im 2Fi{a,b^c, z + 
iO'^) ~ — Im 2Fi{a,b^c, z — iO"*") for z > 1. Furthermore 
Im 2Fi{a,b,c, z) = for z < 1. Thus, integrating sym- 
metrically above and below the real x axis, as per Eq. 74, 
amounts to taking only the real part of 2-^1(0, &, c, z); 
whence we obtain our final result 



2TriiyTai 



.(w) = l + 



p2iujx \ 

/3sinh^)^'^^^ 



2iy/2^ r(i + 2^/i2 

1 

;l-^2^/i2^- 




dx X 



coth ^ ■ 



2' 2 



(75) 



Using the definition of the crossover scale T* = 47r/i^ (see 
Eq. 45), this gives the announced result Eqs. 9, 12 for 
the 2IK model in the special case of perturbation K > 
Kc-, where Saa,a'a' = ^aa'^aa'- By simple extension, 
for K < Kc (corresponding to h < 0) one obtains the 
same crossover but with Saa,a'a' ~ ^S^a'^aa' ■ Finally, 
we note that taking the limit /3 ^ 00 of Eq. 75 yields 
correctly Eq. 53. 

In the next section we show that an identical crossover 
occurs in the 2CK model due to channel anisotropy. 
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V. CROSSOVER IN THE 2CK MODEL 

The NFL fixed point Hamiltonians of the 2IK 
model and the 2CK model have the same basic 
structure. ^*'^^'^^'^*'^^ Although the underlying symme- 
tries of the 2CK model are different from those of the 
2IK model, the free conduction electron Hamiltonian can 
be written in terms of the same MFs in both cases (see 
Eqs. 35-38). The CFT decomposition22 of the 2CK 
model into C/(l) x SU{2)2 x SU{2)2 symmetry sectors 
(corresponding to conserved charge, spin and flavor) , can 
then be expressed in terms of these MFs: the [/(I) the- 
ory with central charge c = 1 consists of a free boson or 
equivalently two MFs x'j (j = 1, 2); the spin SU{2)2 the- 
ory with c = I consists of three MFs Xs = (Xi 7X21X1'") 5 
similarly the flavor SU{2)2 theory consists of three MFs 



X/ 



(X2 ' ~Xi 5 ~X2")- The charge, spin and flavor cur- 



rents can also be written in terms of the MFs correspond- 
ing to those symmetry sectors, as given in Eq. C3 of Ap- 
pendix C. 

In particular, the NFL fixed point Hamiltonian is of the 
form of Eq. 39, with a boundary condition that is again 
simple in terms of the MFs. In the 2CK model, the NFL 
physics arises due to a modification of the boundary con- 
dition in the spin sector only (the free boundary condi- 
tion pertains in charge and flavor sectors). The nontrivial 
boundary condition can be accounted for by defining the 
scattering states Xsix) — ~Xs{~x) and Xf{x) = Xf{—x), 
X'j{x) ~ x'ji^x) for j = 1, 2. Indeed, the finite-size spec- 
trum at the NFL FP^^ can be understood in terms of 
excitations of a free Majorana chain. '^^ 

The NFL fixed point of the 2CK model is destabilized 
by certain symmetry-breaking perturbations. These per- 
turbations can again be matched to MFs, with the cor- 
rection to the NFL fixed point Hamiltonian being of the 
form of Eq. 42 in the simplest case of channel anisotropy 
Ai (X A^ 7^ (see Table I). 22 

Importantly, it was shown recently in Ref. 19 that the 
NFL fixed points of the 2CK and 2IK model are in fact 
identical in the sense that they both lie on the same 
marginal manifold parametrized by potential scattering. 
Indeed, the low-energy effective model for the 2IK model 
in the limit of strong channel asymmetry is the 2CK 
model, ^^'-"^^ but with an additional tt/2 phase shift felt 
by the conduction electrons of one channel. ^^ For con- 
creteness, we consider now a variant of the standard 2IK 
model in which channel asymmetry appears explicitly: 

H21K = Hq + JlSl-sql + JrSr.sqr + KSl-Sr. (76) 

One thus recovers Eq. 2 at the symmetric point Jl = Jr- 
In the limit Jl ^ Jr, the left impurity is Kondo screened 
by the left lead on the single-channel scale T^. At en- 
ergies ~ T^ the right impurity is still essentially free. 
However, it feels a renormalized coupling to its attached 
right lead, and an effective coupling to the remaining 
Fermi liquid baths states of the left lead (which suffer a 
full 7r/2 phase shift due to the first-stage Kondo effect in 



that channel). The relative effective coupling strengths 
between left and right channels can be tuned by the inter- 
impurity coupling K. Tuning K to its critical value Kc 
yields the 2CK critical point;^^ while deviations K 7^ Kc 
correspond to finite channel anisotropy A^ 7^ in the ef- 
fective 2CK model. In consequence, one expects the NFL 
to FL crossover in the two models to be simply related. 
Since the NFL fixed point itself is the same in both 
2CK and 2IK models (up to potential scattering), ^^ and 
because the correction to the fixed point Hamiltonian 
due to the Ai perturbation is the same,22'23 the RG fiow 
along the NFL to FL crossover is identical. To calculate 
the corresponding crossover Green function, we simply 
incorporate the additional it/2 phase shift felt by the left- 
channel conduction electrons into our scattering states 
definition. Using ip'o-hi^) ~ sigii{x)ipa-L, one straightfor- 
wardly obtains 

{i^aaizi)^l,^,{z2))2CK,A.>0 = (V'L (^OV-'Lc' (^2)): 



T,'tzn 



)2IK,K>Ka 



Z2 



2i 



(77) 

in terms of the analytic crossover Green function for 
the 2IK model given in Eq. 60. Following the steps of 
Sec. IV B, the t matrix follows as 



2nii^r' 



l-r^„,[2^zi.r,2^.^J^)-l]. (78) 

This result was obtained at T = in Ref. 26, where the 
correspondence was checked by explicit numerical renor- 
malization group calculation. 

In the 2CK model with A^ > 0, the left lead is more 
strongly coupled, and completely screens the impurity at 
the FL FP on the lowest energy scales; while the right 
lead decouples asymptotically. The physical interpreta- 
tion of Eq. 78 is thus that a Kondo resonance appears 
in the spectral function on the scale of T* in the a = L 
channel, while the resonance is destroyed in the a = R 
channel (hence the dependence on the flavor-space Pauli 
matrix r^„,). 



VI. GENERALIZATION TO ARBITRARY 
COMBINATION OF PERTURBATIONS 

In Sec. IV we considered the finite-temperature 
crossover Green function in the 2IK model due to the 
detuning perturbation K ^ Kc] while in Sec. V we calcu- 
lated the analogous crossover Green function in the 2CK 
model due to channel anisotropy Az ^ 0. In this section 
we generalize the results to an arbitrary combination of 
perturbations in either model. 



A. Flavor rotation in the 2CK model 

Before discussing the full calculation, we motivate the 
general approach by exploiting a bare symmetry of the 
2CK model, in a simple intuitive example. 
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Unlike the 2IK model, the 2CK model possesses a bare 
flavor symmetry (see Eq. 1). The perturbations A^,, Ay 
and Az break this symmetry, but are themselves related 
by rotations in flavor-space. 

A canonical transformation of the conduction electron 
operators of the bare Hamiltonian is defined, viz 



(79) 



such that the unitary matrix U satisfies 
U(A-f]U^ = |A|r^. With the parametrization 

A = |A|(sin6'cos(?!), sin6'sin0, cos6'), one obtains ex- 
plicitly U = exp (|(— sin0r^ -I- cos(/>r'')). It follows 
that 6H2CKiA^,Ay,A,) -^ SH2ck{0,0,A,), with 
A^ = A^ + A^ -f Al (cf Eq. 8). The physical interpre- 
tation is that the impurity couples more strongly to one 
linear combination of channels than the other. Thus 
the perturbations Ai, A2 and A3 in Table I are simply 
related, and their combined effect enters only through 
Xf. In particular this implies only one fitting parameter 
ci for the different components of the perturbation A in 
the 2CK model. 

The physical behavior in the case of arbitrary A/ can 
now be understood in terms of the situation where Ai 
alone acts using the flavor rotation Eq. 79. For example, 

the Green function ((V'fecrL; V'I-'ctl))'^ probes the physical 
channel era with a = L in the original basis. Using the 
transformation Eq. 79 it can be expressed as 



We express the NFL fixed point Hamiltonian in terms 
the free MF scattering states. 






(80) 



in terms of the Green functions in the rotated basis, 
which correspond to those calculated for A^ only, as con- 
sidered in the previous section. It is then easy to show 
that the t matrix for arbitrary Ay is given by 



2TTivTaa,crai^) = ^ + I IT I i'^'^^'^'faa.aai^) - 1 



n) 



in terms of the t matrix T due to Ai given in Eq. 78. 
The simple rescaling of the spectral function discussed in 
Sec. II B and the precise form of Eq. 10 follow immedi- 
ately. 



B. Emergent symmetries 

In this section we make use of the field theoretical de- 
scription of the NFL fixed point for both 2CK and 2IK 
models in the presence of relevant perturbations. '"'^'^^"^'^ 
A large 5*0(8) emergent symmetry at the fixed point al- 
lows these perturbations to be related by a unitary trans- 
formation, in full analogy to the method demonstrated 
explicitly in the previous section for the case of the bare 
flavor symmetry in the 2CK model. 



Hqcp — Hfp[x] + SHqcp, 



(82) 



where Hpp[x] is given in Eq. 39, and with x a- vector of 
the 8 MFs. As already commented, the structure of the 
NFL fixed point Hamiltonian is the same for 2CK and 
21K models; only the definition to the scattering states 
is different. To fix notation, we define 

{Xi,-X8} = {-sign(x)x^,x{,X2:Xi,X2>xf ,Xi,X2} 

(83) 
for the 2IK model; and 

{Xi,---X8} = 

{ -xf , X2 . -x{ , sign(x) Xi , sign(x) X2 , sign(x)xf , Xi , X2 } 

(84) 

for the 2CK model. With this ordering of components, 
the correction to the NFL fixed point due to relevant 
perturbations is given by 



SHqcp oc i^ AjXj(0)a, 



(85) 



with o a local MF as before. The Ai perturbation cor- 
responding to K y^ Kc in the 2IK model or A^ ^ in 
the 2CK model was considered explicitly in Eq. 42. The 
other coupling constants are defined in Table I; with the 
single resulting crossover energy scale being the sum of 
their squares, Eq. 8. For a detailed derivation of Eq. 85, 
see Refs. 17,52,53. The two MFs xTjXs corresponding to 
the real and imaginary parts of the total charge fermion, 
are not in fact allowed in the 2CK and 2IK models due 
to charge conservation. 



C. Unitary transformations 

The crucial observation following from Eq. 85 is that 
only the linear combination A~^ X^^^i "^jXji^) of the 8 
MF scattering states participates in the crossover. The 
particular linear combination depends on the ratios of 
the various perturbations (for example K — Kc, Vlr, Bs 
in the 2IK model). 

From Sees. IV and V, we know the crossover Green 
function caused by the Ai perturbation. The strategy 
is thus to fix the Ai perturbation as the direction in 
the 8-dimensional space of perturbations along which the 
Green function is known, then use an 5*0(8) rotation to 
obtain the general crossover Green function. 

We search now for a unitary operator UU'^ = 1 that 
transforms the full Hamiltonian with an arbitrary com- 
bination of perturbations into one involving the single 
perturbation Ai. Specifically we demand that. 



UHpp[x]U^ ^HfpUI 
USHqcpU^ ^ iXxiiO)a. 



(86) 
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This transformation is accomplished by an operator that 
rotates the S-component vector x in the 8-dimensional 
space of perturbations. The 28 generators of such rota- 
tions are of the form i J dxXjix)Ajj>Xj'ix), where Ajjr 
is a real antisymmetric 8x8 matrix. It is easy to verify 
that the desired operator satisfying Eq. 86 is 



U 



J f^^ dxxiix)xA^) 



(87) 



where 



xA^) ^ ^l^'^>^jxA^)■ 
One can apply this transformation to the 
expectation value of an operator written 
in terms of the original electrons, such as 
the Green function {'^aa{x)'ij?^,^,{x')) hqcp = 

m.o.{x)WU^l^,[x')W)uHc,cpU^-, The^ crucial 

property of the unitary transformation Eq. 87 is that it 
acts as a simple rotation also on the electron fields. This 
occurs due to the existence of linear relations between 
the 28 quadratic forms of the original electron fields and 
of the MFs xf^ '"^ discussed in Appendix C and Ref. 55. 
As a simple relevant example, consider the 2IK model, 
perturbed by a combination oi K — Kc and tunneling 
Vlr, such that only Ai and A2 are finite. In this case 
the unitary operat or reads U = ^^ I^,^ 'i'>'xiix)x2{x) ^-^^^i 
A_L = A2 and A = y^Xf + X^. Using Eq. 83, Xi(a;)x2(a^) = 

— sign(x)x2" (2^)Xi (''c)- The quadratic form X2^(2^)Xi (2^) 
is related to a quadratic form for the original electrons 
X^(a:)x((a;) = —^ip^r^ip, as shown in Appendix C. The 
operator U can now be understood as a simple rotation 
of electron fields, 

C/V'.„(x)[/t = ^ A/i,,,„,V'.'a'(a:), (89) 

a' a' 

for x ^ 0; and where the rotation matrix acts here in 
flavor space. 



(^d,.'.'). 



A/i,.'a')^^,,^^o"^-' 



6aa' cos - T iTaa' Sin -^ 



(90) 



The Green function then follows as 



{i:aa{x)ijl,^,ix'))H. 



QCP / , 



M^a.^,a^(M< )a[a'^^a'a' X (^<Tiai(a;)V'I;„;(a;')>raQc^m, 

(91) 

where x > and x' < is assumed. In terms 
of complex coordinates zi ~ t + ixi and Z2 = 
ix2 (with xi > 0, a;2 < 0), the full Green func- 
tion (VVa(^i)V'o-'a'(-^2))_f/Qcp ^^ then obtained from 



Gaiai,a{a[{^^^f^) = ^ (V'ctiQi (^1 ) V'g./ ^'^ (^2)) ;7Hqcp C/t : ^S 

given in Eq. 60 for the case of finite Ai < in the 2IK 
model. 

Now we define a 4 x 4 unitary Fermi liquid scattering 
S matrix for the 2IK model 

SIZ,^,^-(m>-M<^) (92) 



such that 

{i>c7a{zi)lpl>a'iz2))HQcp = Saa,c7' a'Gaa,. 



2i 



(93) 

Using Eq. 90 one obtains Scra,a'a' = —Saa'{Saa' COS 6 ~ 

i sin 6'r^^,) = ^o-o-/ — ^ °°'^ ' ^"^°°' for the case of finite Ai 
and A2 in the 2IK model; consistent with Eq. 11. 
For arbitrary combination of {Ai, ..., Ag}, one has 



t/ = e 



„i| /^ dx sign(x)^+[A2r="+A3T''-|-(AB-ff)T"liA/Aj_ 



(94) 



From Eq. 89, it then follows that 

'A2r^ + A3T'^ + (AB-a)r^ 



M< = cos - =F * sin - 
2 2 



A, 



(95) 
(suppressing spin and channel indices). Using Eq. 92 we 
recover our final result for the 2IK model, Eq. 11. 

Following the same steps for the 2CK model (and not- 
ing Eq. 77), we obtain 






5^^.^''.' = -(A/>-^M<' ,, (96) 



with 



^ d .. e ( A2r^ - X^T- ± {Xb ■ S)t' 
M< = cos - — I sm - ; 



yielding precisely Eq. 10 for the 2CK model. 



(97) 



VII. NUMERICAL RENORMALIZATION 
GROUP 

Wilson's numerical renormalization group^ (NRG) has 
been firmly established as a powerful technique for the 
accurate solution of a wide range of quantum impu- 
rity problems.^ Its original formulation provided access 
to numerically-exact thermodynamic quantities for the 
Kondo^ and Anderson impurity'"'^ models. An increase in 
available computational resources subsequently allowed 
straightforward extension to multi-impurity and multi- 
channel systems.'^ 

More recently, the identification of a complete basis 
within NRG (the 'Anders-Schiller basis' comprising dis- 
carded states across all iterations'^) has permitted rigor- 
ous extension to calculation of dynamical quantities. In 
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Figure 10: Spectrum ta-L{u), T) vs uj/D for the 2CK model with uj = 0.15 and small channel anisotropy 2z/Az = 10~^ at various 
temperatures T/T* = 10~^, 1, 10, 10^, approaching t^L ~ \ from above. Dotted lines are results from full NRG calculations; 
solid lines are exact results from Eq. 9 for the NFL to FL crossover. 



particular, equilibrium spectral functions can be calcu- 
lated using the full density matrix approach, ^^'^^ yielding 
essentially exact results at zero-temperatures on all en- 
ergy scales. Although discrete NRG data must be broad- 
ened to produce the continuous spectrum, ^^ artifacts pro- 
duced by such a procedure are effectively eliminated by 
averaging over several interleaved calculations (the so- 
called z-trick^''). Indeed, resolution at high-energies can 
be further improved by treating the hybridization term 
exactly.^^ 

Our exact analytic results were tested and confirmed 
by comparison to NRG at T = in Ref. 26. 

Due to the logarithmic discretization of the conduction 
band inherent to NRG,** finite-temperature dynamical in- 
formation cannot however be captured'^^ on the lowest 
energy scales |w| < T. But spectral functions for |w| > T 
are accurately calculated, and the total normalization of 
the spectrum is guaranteed, '^^ implying that the total 
weight contained in the spectrum for |a;| < T can be de- 
duced. From a scaling perspective, one expects RG flow 
to be cut off on the energy scale |a;| = 0{T), so that there 
should be no further crossovers in spectral functions for 
|w| < T. The somewhat arbitrary strategy^^ commonly 
employed is thus to smoothly connect the spectrum cal- 
culated at a; « ztT. in such a way as to preserve the total 
weight. 

Our exact finite-temperature results for the crossover 
t matrix of the 2CK and 2IK models thus offers the per- 
fect opportunity to benchmark NRG calculations for in- 
teracting systems exhibiting a non-trivial temperature- 
dependence of their dynamics. 

For concreteness, we consider now the 2CK model with 



channel-anisotropy A^ > (see Eqs. 1 and 3). To obtain 
the numerical results, we discretize flat conduction bands 
of width 2D logarithmically using A = 5, and retain 8000 
states per iteration in each of z = 6 interleaved NRG 
calculations.^ All model symmetries are exploited. 

To ensure the desired scale separation T* ^ Tk, we 
take representative i/J = 0.15 and small 2zyA2 = 10^^, 
yielding Tk/D = 2 x lO^^ and T*/D = 9 x lO^i'' 
[Tk ^ D exp{^l/i>J) is defined in practice here from 
the t matrix Eq. 6, t„a{u} ~ Tk,T = 0) = ^ T* 
is defined according to Eq. 9, corresponding here to 
taa(w = T*, T = 0) ~ 0.95). From Eq. 8 and Table L we 
thus obtain ci ~ 14. The t matrix for this 2CK model 
can be expressed as. 



' aa^a' a' \^ , -^ } ^aa'^aa' \ ^^ 



Gc,{uj,T) (98) 



with Ja = J ± h^z for a ^ L,R and where 

(99) 
As usual {{A-,B))^rp is the Fourier transform of 

the retarded correlator {{A{ti);B{t2)))T = —iO{ti — 
t2){{A{ti), B{t2)})T- The alternative expression given 
in Ref. 12 is: 



TrV%a,crai(^,T) 



1 



G^{u,T) 



^TTI^JaJ Gc{uJ,T) 

' (100) 
where Ga{uJ,T) = ((^otrc,; V'Jaa>>^,T is the Green func- 
tion for the 'O'-orbital of the a = L,R Wilson chain. ^ 
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Both Ga{oJ,T) and Ga{ui,T) can be obtained directly 
by NRG, but Eq. 100 gives much better numerical 
accuracy, ^^ and is employed in the following. The de- 
sired spectral function t„a[ijj,T) is then obtained from 
Eq. 6, and is plotted in Fig. 10 as the dotted lines for 
temperatures T/T* = 10"\ 1, 10, 10^, as in Fig. 3. The 
corresponding exact results for the NFL to FL crossover 
from Eq. 9 are plotted as the solid lines. As immediately 
seen, near-perfect agreement is obtained for all energies 
\(jj\ <?C Tr- and temperatures T <^ Tr- where comparison 
between numerical and exact results can be made. 

To obtain such an agreement, we found that high- 
accuracy NRG calculations must be performed. In par- 
ticular, the region |a;| ~ T was most problematical, with 
artifacts only being removed upon averaging over several 
band discretizations, and necessitating a large number 
of states to be kept at each NRG iteration. The precise 
shape of the numerically-obtained spectrum then still de- 
pends on how the discrete data is smoothed. We found 
that the broadening scheme described in Ref. 59 pro- 
duced the best results: for z = 6 and A = 5 as used here, 
a broadening parameter b = 0.25 and kernel-crossover 
scale luq = T/1.5 were optimal. 

It should also be noted that if the correction factor'^ 
Ajy = 5^3Ylog(A) is used in the NRG calculations (in 
which case J^ — > JqAa), then the many-particle energies 
used to calculate the density matrix must be accordingly 
scaled (i?Ar(r) — >■ EN{r)A\) so that the results are in- 
dependent of the discretization parameter A and hence 
approximate accurately the desired A = 1 limit. 



VIII. OTHER EXACT CROSSOVER 
FUNCTIONS 

As discussed in the previous sections, the Hamilto- 
nian controlling the NFL-FL crossover in the 2CK or 
2IK models has a free fermion structure in terms of MFs. 
In fact, this feature allows calculation of various quan- 
tities along the crossover. The difficulty of such calcu- 
lations is dictated by the relation between the physical 
quantity of interest to the MFs. In the preceding sec- 
tions we concentrated on the two-point function of the 
electron field (the Green function), related^^ here to the 
one-point function of the magnetization operator in the 
Ising model (which is in turn related non-locally to the 
Ising MFs^^). More generally, 2p-point functions of the 
electron field are related^^ to p-point functions of the 
magnetization operator. Multi-electron correlators can 
thus in principle be calculated in this way, but require 
knowledge of the corresponding multi-point correlation 
functions of the Ising magnetization operator. 



A. Impurity entropy 

Since thermodynamic quantities are local in the MFs, 
their calculation is rather straightforward. Here we will 



focus on the NFL-FL crossover of the impurity en- 
tropy, following closely the earlier calculations for the 
2CK modeP"^'^® and the 2IK model, ^"^ performed in the 
Toulouse limit. The Toulouse limit corresponds here to 
maximal spin anisotropy in the exchange couplings, and 
as such breaks the overall SU{2) spin symmetry of the 
models. Although the high-energy ^ Tk crossover to 
the NFL fixed point is strongly affected by large spin- 
anisotropy, we stress that for low energies <^ Tk (and 
given a clear scale separation T* ^ Tx), the results be- 
come formally exact, and are universally applicable to 
the SU{2) symmetric case of interest. 

The key point is that the spin-anisotropy perturbation 
is RG irrelevant at the NFL fixed point. In particu- 
lar, the effective theory obtained in the Toulouse limit 
describing the NFL-FL crossover due to relevant pertur- 
bations such as channel anisotropy or magnetic field in 
the 2CK model^^ or detuning K — Kc, staggered mag- 
netic field, or left-right tunneling in the 2IK model, ^^ act 
exactly as in Eq. 85. A detailed discussion for the 2IK 
model can be found in Ref. 52. 

Turning now to the crossover in the impurity entropy, 
one finds'^^ that 



5(T) = ilog(2)-K5 



in terms of the universal function 



(101) 



Sit) = i 



^<^7 



1 



log 



1 



, 1 1 
.^2 + 7 



(102) 
defined in Ref. 28 for the limit Tk -^ oo. Here, ipi^) 
is the psi (digamma) function and T* is a particular 
definition of the NFL-FL crossover scale (proportional 
to our definition, Eq. 9, such that T* ~ y x T* with 
y « 4.6). Two regimes can thus be distinguished. In 
the FL regime, obtained for T <ti T* , the impurity is al- 
ways completely screened: S ^ -^ i^j. By contrast, 

in the NFL regime, T* <^ T <^ Tk, the impurity en- 
tropy is close to ilog(2). Interestingly, we find that in- 
dependently of the relevant perturbations which act, the 
entropy crossover is always given by the universal func- 
tion Eq. 101 in the hmit T* < Tk, in both 2CK and 
2IK models. This is of course not the case for the Green 
function, because the FL scattering S matrix is affected 
differently by different perturbations (see Eqs. 9-12 and 
Figs. 3, 5). 

The 2CK model has also been solved exactly using the 
Bethe ansatz,^ yielding the full evolution of thermody- 
namics in any parameter regime. However, it cannot be 
seen directly from the Bethe ansatz equations that there 
is an emergent SO (8) symmetry at the NFL fixed point, 
or that this leads to a single NFL-FL crossover function 
for the entropy, Eq. 101, regardless of the perturbation 
causing the crossover. Indeed, the fact that the same 
crossover occurs in the 2IK model cannot be extracted 
using Bethe ansatz since the 2IK model is not integrable. 
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Figure 11: Impurity contribution to entropy S'imp(T) vs T/T* for the 2CK model (upper panels) and the 2IK model (lower 
panels) in the presence of various perturbations. Entire temperature-dependence calculated by NRG for the full models (black 
lines); low-temperature T <^ Tk behavior in each case compared with the single exact NFL-FL crossover function Eq. 101 
(red circles). All results presented for uj — 0.25. Left: effect of channel asymmetry Aul^^ = ±10^''' (2CK) or deviation from 
critical coupling {Kc - K)/D = ±10"^ (2IK). Center: effect of including also finite left/right tunneling, 2v/^^ = ±10~^ (2CK) 
or 2vVlr = ±10"'' (2IK), with 4i^A2 = {Kc — K)/D = ±10"^ as before. Right: effect of including finite magnetic field, 
B'/D = ilQ-^/^ (2CK) or Bl/D = ±10"^/^ (2IK), again with 4vA^ = (A'c - K)/D = ±10"^ Parameters chosen to allow 
direct comparison to Fig. 1 of Ref. 26. 



In Fig. 11 we present numerically-exact NRG results 
for the temperature-dependence of the entropy due to 
various perturbations in the 2CK and 2IK models to con- 
firm the validity of the field theoretic description. 

As in Fig. 10, we exploit all model symmetries to ob- 
tain high-quality numerics, discretizing fiat conduction 
bands of width 2D logarithmically, using A = 3 here, 
and retaining 8000 states per iteration in a single NRG 
calculation.^ 

At low temperatures T ^ Tk (and since T* ^ Tk), 
we obtain an essentially perfect agreement between the 
exact result Eq. 101 (points) and NRG data (solid line). 



B. Non-equilibrium transport in two lead devices 



tichannel setups have been considered in Refs. 64,65. 

We consider a finite source-drain voltage V across left 
and right metallic leads, which are exchange-coupled to 
a single impurity spin. To this system we add a small 
but finite channel anisotropy perturbatioir, correspond- 
ing left/right tunneling mediated via the impurity spin. 
The setup is illustrated in Fig. 12. The corresponding 
Hamiltonian is given by Eqs. 1 and 3, with finite A^; and 
possibly magnetic field B, but irow with left/right lead 
chemical potentials at ±V/2. 

The applicability of our exact solution is in the paranr- 
eter regime A^ ^ J, so that the system is close to the 
NFL critical point. This situation is not in practice ob- 
tained in standard quantum dot devices, although more 
sophisticated experimental techniques such as those em- 



It was shown in Refs. 52,53,63 that the effective free 
fermion theory of the 2IK model allows to calculate cer- 
tain non-equilibrium quantities. Finite conductairce was 
found to arise iir the weak coupliirg linrit of 2IK systems 
close to the critical point T* <C Tk at low energies <C T* . 
This result was understood in terms of the growth under 
RG of the left-right tunneling perturbation Vlr- Here, 
we generalize these results to the 2CK model, which has 
the sanre effective free fernrion description. Related nrul- 




Figure 12: Schematic illustration of a non-equilibrium 2CK 
setup. We consider the case of A^: <^ J. 
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Figure 13: Scaling function for the nonlinear conductance of a 
2CK device with small A^,. From top to bottom at the peak: 
T/T' =0,0.25,0.5,1. 



ployed in Ref. 24, do allow suppression of cotunneling 
perturbations such as A^,. 

As per Eq. 8, the crossover energy scale is T* = 
{cii^Aj^fTx + \cbB\^/Tk- In the limit where vA^ 
is initially very small, we thus have T* ^ Tr-. At 
higher energies > Tk, we then expect conductance to 
be very small ex {vAx)'^, corresponding to the weak cou- 
pling limit. However, upon reducing the energy scale 
E = max{ey,r} below Tr-, the conductance starts to 
increase since A^, switches on a relevant operator with 
scaling dimension 1/2 near the NFL FP. Below T* , a 
characteristic peak in the conductance is thus expected, 
signaling growth of the relevant operator to order one. 

The exact lineshape of the non-equilibrium conduc- 
tance peak can be calculated from the fixed point Hamil- 
tonian, Eq. 82 (including the correction due to relevant 
perturbations given by Eq. 85). The dependence on the 
ratio between the magnetic field and the tunneling per- 
turbations is obtained using the 5*0(8) rotation outlined 
in Sec. VI. The method of calculation, and the result for 
the universal lineshape of the peak, was obtained for the 
spin-exchange anisotropic version of the model by Schiller 
and Hershfield in Ref. 66. As argued in the previous sub- 
section, we can borrow those Toulouse-limit results if we 
restrict attention to the low energy crossover. The final 
result for the nonlinear conductance is thus 



G^GqF 



T eV 



, , 1 /I 1 iv 

F i, v\ = Re Vi - + 1 

^ ' ^ 47ri ^ I 2 47ri 27ri 
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with Go = ^1 — *p~^ ^^^ '/'I ^^^ trigamma function. Note 
that the definition of T* here is as in Ref. 52. At T = 0, 
one obtains G/Gq = [l + {2eV/T*f]-'^] while at zero-bias 
V ~ and low-temperatures T ^ T*, the asymptotic 
conductance is G/Gq ^ 1 - {2TrT/V3T*)^. The fuU bias- 
dependence of conductance for various temperatures is 
shown in Fig. 13. 



IX. CONCLUSIONS 

In this paper we present a rare example of an exact 
nonperturbative result for the finite-temperature dynam- 
ics of a strongly correlated quantum many-body system. 
We focus on the two-channel Kondo and two-impurity 
Kondo models; although the same low-energy physics 
characterizes a wide class of quantum impurity problems 
in which competition between two conduction channels 
causes a frustration of screening. The unusual non-Fermi 
liquid critical points of these systems are destabilized 
by various symmetry- breaking perturbations, naturally 
present in experiment. In consequence, a crossover to 
regular Fermi liquid behavior always occurs on the lowest 
energy scales. Exploiting the connection^'^ to an exactly- 
solved classical boundary Ising model, •^"^■'^^ we calculated 
the exact finite-temperature crossover Green function. In 
quantum dot systems which could access this crossover, 
the relevant experimental quantity is conductance, which 
we extract from the exact Green function. 

Remarkably, we show that due to the free fermion 
structure of the effective low-energy theory in terms of 
Majorana fermions and a large emergent SO{8) symme- 
try, a single universal function pertains for any com- 
bination of perturbations in either model. This single 
crossover is also starkly manifest in the behavior of ther- 
modynamic quantities such as entropy; as confirmed di- 
rectly by NRG. 

The method developed in this paper goes beyond the 
impurity models we considered explicitly, and finds pow- 
erful application to a wider family of systems. At heart, 
our solution relies upon a formal separation of the the- 
ory into a sector containing all the universal crossover 
physics, and a sector acting as a spectator along this 
crossover. Importantly, the crossover is confined to a 
sector which can be identified with Ising degrees of free- 
dom, described by a minimal conformal field theory with 
central charge c = 1/2. For example, in the two-channel 
Kondo models studied here, the full set of degrees of free- 
dom consist of a c = 4 GET, but a large c = 7/2 sector of 
the theory plays no role in the crossover from non-Fermi 
liquid to Fermi liquid physics. 

Interestingly there exist other models (whose full set 
of degrees of freedom are not necessarily described by a 
c ~ 4 GET) which undergo precisely the same crossover 
due to their underlying c = 1/2 Ising sector. Those in- 
clude certain Luttinger liquids containing an impurity,'^ 
and coupled bulk and edge states in certain non-Abelian 
fractional quantum Hall states^"^'^® (see also Ref. 69). 

There are further interesting generalizations and ques- 
tions arising from this work. For example, the two- 
channel Kondo effect evolves continuously as interactions 
are switched on in the leads, as was shown in the case 
of Luttinger liquid^°'^^ and helical liquid^^"^'* leads. It 
is an open question as to whether the low-temperature 
crossovers in the presence of such interacting leads are de- 
scribed by the same boundary Ising model, or e.g. by cou- 
pled boundary Ising models. It would also be interesting 
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to use the present formulation of the crossover in terms 
of a minimal Ising theory to study time-dependent phe- 
nomena, quench dynamics, and other non-equilibrium 
physics. 
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Figure 14: Diagram for the inelastic contribution to the t 
matrix. It describes interaction between two fermionic species 
A ^ B. In the ICK model A,Bg {t,i}; in two-channel 
Kondo models A,B e {t L,i L,^ R,i R}. 



Appendix A: Perturbation theory around the FL 
fixed point 

In this appendix we consider the 2IK model for K < 
Kc, and its FL fixed point describing the ground state 
where each impurity forms a Kondo singlet with its at- 
tached lead. In particular, we calculate the t matrix for 
uj,T ^ T* as a stringent consistency check of our full 
crossover t matrix, Eq. 9. Indeed, we also see that the 
multiplicative function f{(3h^) that we included in Eq. 54 
is precisely needed to reproduce the correct FL limit. 

We use here the Fermi liquid theory of Nozieres,^^ ap- 
plied to the 2IK model. Naturally, our derivation of the t 
matrix in the vicinity of the FL fixed point follows closely 
the analogous calculation for the simpler single channel 
Kondo (ICK) model. Thus, we first recap some of the 
basic concepts and results for the ICK model. •^^''^^'■^* 

The irrelevant operator in the effective Nozieres 
Hamiltonian'*^ for the FL fixed point of the ICK problem 
may be written in CFT language as^^ 



SH 



ICK 






(Al) 



2T, 



where J{x) = ip^]y{x)^-^il'' „,{x) is the spin current for a 
single channel of conduction electrons (implicit summa- 
tion over repeated indices is implied). The first term of 
the second line may be interpreted as an elastic single 
particle scattering, and the second term can be inter- 
preted as a residual electron-electron interaction giving 
rise to inelastic scattering. Accordingly, one can sepa- 
rate the contributions to the t matrix as Tick iw,t) = 
Tel {w, t) + %n {w, t), where22.54 



-TTl/Tel {w, t) = i — w — iw'^, 



2,2i 



-TTvTiniWjt) = --[w +TT t 



(A2) 



[For simplicity we omit spin indices; TicKaa' (w, t) = 
&aaTicK {w,t)\. Here w = uj/T'j^, t == T/T'j^, and T^ is 
a particular definition of the Kondo temperature.^^ The 
fermionic diagram yielding the inelastic contribution is 
shown in Fig. 14. 



The imaginary part of the t matrix in the ICK model 
can thus be expanded as 



Im Tick {w, t) = Im Tick (0, 0) + aW 
where from Eq. A2 one obtains 

a/fo|iCK = 3/7r2~ 0.3039. 



bt' 



(A3) 



(A4) 



Similar to the Nozieres Fermi liquid theory for the 
ICK model, an effective Fermi liquid Hamiltonian was 
constructed using CFT methods based on the emergent 
S0(7) symmetry of the crossover for the 2IK model. ^'^■^^ 
The leading irrelevant operator takes the form^"^ 



/t 
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Jfj — 6JlJb.)x=o 



(A5) 



where Ja{x) = ilj^l^^{x)-^tp'^^,{x), ip' is a scattering 
state incorporating the 7r/2 Kondo phase shift, ip' aai^) — 
sign(x)ipaa{x) ^ and T* is a particular definition of the low 
energy crossover scale. ^'^ We now calculate the t matrix 
resulting from this Hamiltonian in the FL regime w, T <C 
T* . In the following we suppress the indices cr, a, and 
note that the t matrix is proportional to d^a'^aa' hi the 
present situation. 

Comparison of the irrelevant operators in the ICK 
and 2IK models, Eqs. Al and A5, shows that the first 
two terms in Eq. A5 are identical to the Nozieres ir- 
relevant operator (up to the exchange of energy scales 
T* <-)• Tk)- In consequence the elastic and inelastic 
scattering contributions within each channel of the 2IK 
model are the same as those arising in the ICK model. 
Indeed, from the diagram Fig. 14 we see that to sec- 
ond order they do not yield any mixed terms. Thus 
T21K (w, t) = Tick {w, t) + Tlr {w, t), where Tlr {w, t) 
originates from the third term in Eq. A5, representing 
interaction between channels. We separate the latter into 

Jl ■ Jr = J2a=x y z'^lJr ^^'^ ^^te that the second order 
a — x^y, z contributions are equal, since the quantity of 
interest is invariant under spin rotations, and the Hamil- 
tonian is also SU{2) spin-symmetric. The latter can be 
written in terms of fermion fields as 






2T 



-(V''Ii^LV''.L)(^'Ii^^a>'.fl)- 



(A6) 
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Considering now the t matrix for a single electronic 
species with quantum numbers A = a, a, Eq. A6 de- 
scribes the interaction with a second species of either B = 
a, a or a, a, and thus contributes a term proportional to 
the inelastic contribution %„ iw,t). In fact, the ampli- 
tude for this interaction is zt^^nr, identical in absolute 
value to the intra-lead interaction amplitude between up 
and down electrons [the second term of the second line in 
Eq. Al). But since the contributions to inelastic scatter- 
ing are of second order (see Fig. 14), the sign of the scat- 
tering amplitude is unimportant in calculation of the t 
matrix. Summing over the second species B yields an ex- 
tra factor of 2, yielding Tlr {w, t) = 3T£fi = 6Tin {w, t). 
Putting all the contributions together, we have 



T2IK (if, t) ^ Tel (W, t) + 7Tin (w, t) 



(A7) 



As a result, one again obtains Im T21K {w,t) 
Im T21K (0,0) + aw^ + bt^, but with 



a/b\ 



2IK 



(1 + 2/7)/7r2 ~ 0.13027. 



(A8) 



This result is in perfect agreement with a numerical eval- 
uation of our full finite-temperature crossover t matrix, 
Eq. 9 in the hmit T,uj <$:T*; as demonstrated in Fig. 4. 
This calculation also confirms the need for the function 
/(/3/i2) used in Eq. 54. 

The asymptotic result Eq. A8 also follows from renor- 
malized perturbation theory calculations presented re- 
cently in Ref. 76 for a related 2IK model. 



Appendix B: Perturbation theory around the NFL 
fixed point 

1. General structure 

In Sec. HID we used the analyticity of the Green 
function and the locality of the t matrix to argue that 
Gaa.a'a'i^i ~ ^2) dcpcuds Only ou the difference (zi — 22)- 
This result should hold to all orders in perturbation the- 
ory, as shown explicitly in this appendix. 

Our starting point here is the NFL FP Green function 
for the 2IK model, Eq. 41, written as. 



Gaa.a'a'i^l-i^-'i) —5crcr'5. 



2ir R 



au ^aa 



^sin[^(zi-Z2)]y 

x((Tl(zi)(Tl(z2)), 



(Bl) 



whose factorized form originates from the Bose-Ising de- 
composition of the 2IK model into spin, isospin and Ising 
symmetry sectors. ■^'^ 7 of the 8 MFs (those correspond- 
ing to the part of the fermion field carrying spin and 
isospin quantum numbers) remain free at the NFL FP, 
and thus give rise to the first factor (see Eq. 28). The sec- 
ond correlator involves the chiral part o'i(z) of the Ising 
magnetization operator ai^z) = aL(z)(TL{z*), due to the 
remaining MF in the Ising sector (see Fig. 9). 



When the detuning perturbation Ai oc {K — Kc) acts, 
the NFL FP is destabilized. The perturbation appears 
as a correction to the action, ^^ 



SS=^ Xl f dT€{0,T). 



(B2) 



where e is the CFT d = 1/2 boundary operator from 
the Ising sector interpreted as a boundary magnetic 
field. ^'^ The resulting corrections to the Green function 
can then be calculated within perturbation theory. The 
full crossover Green function is then G^^^;^f (zi, Z2) = 
^aa^'a'i^i ~ ^2) + J2n=i ^nGccm' a' {zi, Z2) , whcrc the 
A'^-th order correction is given by (suppressing spin and 
channel indices) 



5nG{zi,z2) oc Al 



N 






N 



{aL{zi)aL{z2)X{e{Q,T,)). 



(B3) 



Generically, correlation functions up to 3-point func- 
tions are determined by CFT. However, the Ising GET 
is special because essentially all correlation functions are 
known exactly. In particular, Ardonne and Sierra ob- 
tained explicit expressions^* for the correlators appearing 
in Eq. B3. In the case of even N their result reads, ^® 



N 



{cTLizi)aLiz2) JJe(0,Tj)) ex (zi - Z2)" 



1/8, 



]=i 






Hf. 



ijei- 



{Zl - Z2)l^^l 



(B4) 



where the sum is over all subsets of {1, 2, ..., A^}, con- 
taining an even number of elements |/|. Ilf(Af) denotes 
the Haffnian of a symmetric N x N matrix, and is given 

by Hf(A/) == 2"/^(\f/2)! S'^gSiV U^Jl ^U2^-l),a{2i)' ^^^^ 

a a permutation. The set / (containing |/| elements) 
is equal to {1,2,...,A^} \ /. Using the conformal map- 
ping Eq. (59) from the plane to the cylinder, each co- 
ordinate difference z — z' in Eq. (B4) is replaced by 



— sin 



^iz-z') 



The dependence on zi and Z2 is 

through factors which explicitly depend on zi — Z2, and 
terms inside the square-root of the form 



1 



rije/Sin (^(zi - r^-)) sin (|(z2 - Tj) 
With the aid of the trigonometric identity, 
2sin(^(zi-r))sin(|(z2-T)) 
= cos(^|(zi -Z2)j -cos(^|(zi + Z2-2t) 



(B5) 



(B6) 
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and by shifting all tj integration variables by (zi + Z2)/2 
into the complex plane (which can be done without en- 
countering any singularities), the resulting integral in 
Eq. B3 then depends only on zi — Z2. The same con- 
clusion is reached for odd-A^ by a similar calculation. 

Thus SNGcra,a'a'izi,Z2) = 5nG ca,a' a' {zi ~ Z2) for aU 

N , and hence the full crossover Green function depends 
only on [zi — 22). Analyticity of the Green function 

{i'aa{x)w'^,„,{~x)) = Gaa.a'a'ix) thus allows determi- 
nation of {■llJaa{zi)lpl,^,{z2)) = Gcra^a'a' (^^^) "l termS 

of general coordinates Zi and Z2 by analytic continuation. 



2. Leading order perturbation theory 

Here we derive the NFL coefficients /3', S and /?" of 
the asymptotic t matrix discussed in Sec. II B. Since the 
Green function vanishes at the NFL fixed point, the lead- 
ing correction arises to first order. We now use the first 
order result for (cr(a;)) derived in Ref. 43, 



(a(x))W ^h. 



47r 
13 



1/8 



sinh^ 




Oih^). 



2K\; 



Note that 2^1 (51 5; l^^) = ^^7^ with K the complete 
elliptic integral of the first kind. We will also use the 
short distance a; — > limit of this formula. 

{a{x))h,p = -2"/8/ix3/8[log(a;) + 0(1)], (B8) 

vahd for x <^ j3, h^^. Eqs. 44 and B8 give 

h 1 



G{x ^ 0) = - 



7rV2i 



i Vx 



= [log(x)+0(l)]. (B9) 



To obtain the expansion of the t matrix at T = at large 
uj we use G{x -> 0) in Eq. 74. Recalling that T* = Anh'^, 
and writing logx = log(ajx) — logw and y ~ ljx, Eq. 74 
becomes 



t= -+ ^ .„„ 



2,,ylOgy_]OgW 

2 (2^)3/2 ^"V^^-^"^ ^ ■ 



■Im 



dy e 



(BIO) 



Note that the y integral should be made symmetrically 
around the branch cut, as described in Sec. IV B. One 
then obtains Eq.l7a with 



P' = - 
5 ^ - 



14- TT 
\/2^3/2 7o 
I 

72^' JO 



dy 
dy 



s\TL{2y)\ogy 
sin(2;/) 



(Bll) 



with the numerical values of these integrals announced 
below Eq. 17a. At finite temperature and a; = we use 



Eq. B7 in Eq. 44. Taking the hmit w ^ of Eq. 74 and 
defining y = x/ P, we obtain Eq.l7b with 



/?' 



„ _ 4V2 [°^ ^^ y^ (- 



dy 



— coth(27ry) 
2 



sinh(27rj/) 



(B12) 



Appendix C: Linear relations betw^een quadratic 
forms for original fermions and Majorana fermions 

There are 8 x 7 = 28 independent quadratic forms in- 
volving f/'J-Q, and V'ct'q', which together comprise the gen- 
erators of the SO{8) symmetry group. These generators 
are linearly related to the 28 quadratic forms of the MFs 
xf- In this appendix we gather and re-derive some of 
these relations, which can also be found in e.g. Ref. 47. 

First we define a convention relating the Klein factors 
Fa for the new fermions ipA, to the Klein factors Fa-a 
for the original fermions tpa-a- The relations are fully 
determined by"*® 



fIfI = fI F, r , FxFl = fLf, 



X^ s 



tL^lL 



t-R-f^^-R' 



(CI) 
(B7) and by the anticommutation relations'** {Fa,Fb} = 
26 AB, FaF{ = f\Fa = 1. 

We now consider the instructive example of the op- 
erator iXiXi : and use Eqs. 37 and 38 to relate it to a 
quadratic term involving the original fermions: 



F].F\e"^ !+'*'' + LFfF],e''''''f+"'"' +R.c 



2 f "" 



(C2) 



Using Eqs. CI and 36 we obtain. 



*x{xf = --Fl^F^ne'^r--''f'r---Fl^F^Le''^^--''^^-+li.c. 
Finally using the bosonization formula Eq. 35 we have 



ixixi = -o^'tiV'Tfl - oi'UAL + H.c. = Tp''^—ip 



2^^ 



9n-R'' 



In a similar fashion, all of the relations between quadratic 
forms can be determined. Conserved currents in the 2CK 
and 2IK models can be expressed in terms of the original 
fermions or the MFs, and the relations between them are 
needed for our generalization of the crossover t matrix 
to arbitrary perturbation, as considered in Sec. VI. The 
conserved currents of the 2CK model are, 

charge : J = -■i/'V = «X2Xi 

spin : X =- 2^^^^' "" ~*(X2xf :Xf Xi,X?X2) 
flavor : Jf = -V'V-0 = (-Jx{xf ,«xf X2'«X2x{)- 
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Equivalently. one can define a 3-component spin vector 

Xs = ixhxh Xi) and flavor vector Xf = {x^, -x{, -X2) 
such that, 

— > — 2 — ' — Z 

Js = -^Xs X Xs, J/ = —X! X XS- (C3) 

Furthermore, the 9 spin-flavor current components can 
be expressed as, 



1 



^t^VV' = *(x.r(x/)', (a,6 = x,y,z). (C4) 



Thus, the decomposition of the 2CK model into U(\) x 
SU{2)2 X SU{2)2 charge, spin and flavor sectors can be 
understood also in terms of MFs. 

In the 2IK model, there is no flavor SU{2) symmetry 
since each channel couples to a different impurity. How- 
ever, one can make use of the SU{2) total spin current 
Js, as well as SU{2) isospin currents for each channel II 
and Ir, where 

^" " 9 XI ^itaV'.ra, I^ = V'taV^^c- (C5) 



In terms of MFs, we have 



JZ _ JZ _ J 



z _ JZ _ -J^J 



f =«X2Xi, 



and 



If+I% 



f. 



+ ^]? = «xix2 



tX tX 



ixixi- 



(C6) 



(C7) 



Hence one can understand the conformal embedding of 
the 2IK model as a SU(2)2 x SU{2)i x SU{2)i x Ising 
decomposition into total spin, left/right channel isospin 
and Ising sectors. 3 of the 8 MFs represent the total 
spin sector; 4 represent the isospin symmetry sectors (the 
charge and flavor MFs); and the remaining MF, X2 ' i^ 
associated with the Ising model (and restores the total 
central charge c = 4) . 
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